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Spectral parameter (equation 6) 

Magnetic field strength (equations 1, 9) 

Energy (equation 1) 

Energy limit (equation 10) 

Energy threshold (equation 12) 

Pitch angle dependence functions (equations 6, 7) 
Average storm flux (Figure 64) 

Storm model parameters (Figure 64) 

Unidirectional electron flux (equation 1) 

Average differential flux in channel i (equation 10) 
Omnidirectional flux (equation 12) 

Electron content parameter (equation 8) 

Mcllwain shell parameter (equation 1) 

Pitch angle parameters (equations 6, 7) 

Ratio of Starfish to total flux (equation 23) 
Probability of storm flux (Figure 64) 

Storm ratio (equation 11) 

Solar cycle ratio (equation 15) 

Storm parameters (Figure 64) 

Starfish flux cutoff time (equation 21) 

Universal time from solar minimum (equation 1) 
Spectral hardness parameter (equation 4) 

Local pitch angle (equation 1) 

Equatorial pitch angle (equations 2 y 6) 

Starfish pitch angle cutoff (equation 22) 

Starfish decay time (equation 20) 

Pitch angle cutoff (equation 6) 
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1. INTRODUCTION 


This report is a description of the work performed in the develop- 
ment of the inner radiation zone electron model AE-5. It will be dis- 
tributed to the users of this model together with a document providing 
a brief user-oriented description of the model and associated computer 
programs (Teague et al., 1972). 

A complete description of the omnidirectional flux model is given 
for energy thresholds E 7 in the range 4.0 > Ej/ (MeV) > 0.04 and for L 
values in the range 2.8 > L > 1.2 for an epoch of October 1967. Con- 
fidence codes for certain regions of B-L space and certain energies 
are given based on data coverage and the assumptions made in the anal- 
ysis. The electron model programs that can be supplied to a user are 
referred to. One of these, a program for accessing the model flux at 
arbitrary points in B-L space and arbitrary energies, includes the 
latest outer zone electron model and proton model. Efforts have been 
made to ensure that the two electron models are compatible at the inter- 
face of L = 2.6. In addition, a program for integrating the model flux 
along an arbitrary orbit can be supplied. 

The model AE-5 is based on data from five satellites, 0G0 1, OGO 3, 
1963-38C, OV3-3, and Explorer 26, spanning the period December 1964 to 
December 1967. Because these data primarily provide unidirectional 
fluxes, a unidirectional flux form of AE-5 was first derived. This 
unidirectional form can be supplied to users in a variety of ways, as 
described in this document. The method used to derive AE-5 is dependent 
on both L and energy, owing to the varying effect of magnetic storms in 
the inner radiation zone. For most L values at energies below 700 keV, 
a sufficient average description of the inner zone flux can be obtained 
without regard to storm effects, and for this region an analytical de- 
scription of the flux is presented. The derivation of this analytical 


1 



form is described in Sections 2 and 3. At high L values and high 
energies , however, the effect of storms on the average flux becomes 
significant and is included in AE-S. The derivation of the high-energy 
portion of the model is discussed in Sections 4 and 5. In the course 
of deriving the analytical quiet day (i.e., storm free) flux, two 
temporal variations were modeled. First, long-term increases in the 
inner radiation belt flux were observed that were considered to be 
correlated with increased solar activity over the period 1964 to 1967. 

A tabular form of this solar cycle effect is presented for energies 
E < 700 keV. At higher energies similar temporal variations are ob- 
served, but these generally become insignificant in the context of an 
average flux because of the increased significance of magnetic storms. 
Solar cycle effects are discussed in Section 6A. The second temporal 
variation modeled was that caused by the decay of residual electrons 
from the Starfish nuclear explosion of July 1962. A model of these 
residual electrons is presented that describes the decay time as a 
function of energy and L value. In addition, cutoff times are pre- 
sented that represent the times at which the flux has essentially 
decayed to the quiet day background. For the AE-5 model with epoch 
October 1967, Starfish electrons remain significant only at low L 
values and intermediate energies (E = 1 MeV, approximately) . Starfish 
electrons are discussed in Section 6B. 
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2. THE ANALYTICAL UNIDIRECTIONAL FLUX 

The differential unidirectional electron flux in the inner radia- 
tion belt is a function of five variables: the local magnetic field 

strength B, the local pitch angle a, the local L value, the particle 
energy E, and the universal time T. That is, 

j = j(a,B,E,L,T) (1) 

The first adiabatic invariant relates the local pitch angle and field 
strength to the equatorial pitch angle a 0 and the equatorial field 
strength Bq. The equatorial field strength is a function of L and 
the magnetic moment of the earth only. Equation 1 then becomes 

j = j(a 0 ,E,L,T) (2) 

which is referred to as the equatorial pitch angle distribution. An 
alternate form of equation 1 is the perpendicular flux. The local 
pitch angle and field strength can be reduced to the mirror point 
field strength B x by using the first adiabatic invariant, and equation 
1 becomes 


j = j(B lf E,L f T) C3) 

A complete description of the unidirectional flux at given E, L, and T 
is then given by the equatorial pitch angle distribution for 90° > a Q > 
a c or by the perpendicular flux distribution for B c > B x > Bq, where 
a c and B c are, respectively, the equatorial pitch angle and the field 
strength corresponding to the atmospheric cutoff. 
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Both equatorial pitch angle distribution and perpendicular flux 
data sets are used in the development of the present model 
(Section 3). The analytical model is given in terms of a spectrum, an 
equatorial pitch angle dependence, and sets of coefficients describing 
the temporal changes occurring over the time period September 1964 to 
December 1967. It is considered that the total electron flux in the 
inner belt at any given time is composed of four components: (1) a 

quiet day (that is, periods remote from the effects of magnetic storms) 
background flux at solar minimum, (2) a quiet day component that is de- 
pendent on solar cycle, (3) storm time flux, and (4) residual flux from 
the 1962 Starfish nuclear explosion. Since the diffusion equation is 
linear and since the particles already at a given L value will not affect 
the particles diffusing to that L value to a first approximation, these 
components can be separated and studied as separate processes. The com- 
ponents are not necessarily attributed to different physical phenomena 
but are merely modeled separately. It is evident, for example, that the 
flux increase with solar cycle is physically associated with the storm 
time component. The last two components, (3) and (4), are discussed in 
detail in later sections (Sections 6A and 6B) , because the present ana- 
lytical model is primarily an attempt to describe the first two compo- 
nents. Evidently the four components of the flux are not easily sepa- 
rated, and any possible separation is heavily dependent on the L value, 
the particle energy, and the universal time. For instance, in late 
1964 for L < 2, the dominant flux component is the residual Starfish 
flux, and no description of the other components can be obtained directly. 
In late 1967, however, the residual Starfish flux for E < 1 MeV has de- 
cayed to an insignificant level for most L values, and observations of 
the remaining three components can be made. The present analytical 
model describes the inner belt quiet day flux for energies E < 690 keV 
at an arbitrarily chosen epoch of October 1967. No attempt is made to 
separate the first two flux components because no data are available at 
solar minimum, although an elementary description of the solar cycle 
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effect is given for certain L values in Section 6A. At energies 
E > 690 keV, the data available were limited and an analytical quiet 
day model is not presented. Detailed discussion of high-energy elec- 
trons is presented in Section 4. 

The assumed spectral form is that previously discussed by Teague 
and Vette (1971) , and the unidirectional flux at some reference time 
T r , reference equatorial pitch angle a r , and L value is 

j r = a r (L,ot r ,T r ) E exp [-E/X^L'd^T,)] (4) 

where a^. and X r are the reference parameters. Teague and Vette (1971) 
showed that this spectral form provided a better description of the 
quiet day flux than a power law, an times exponential form, or the 
more normally assumed simple exponential law. The reference time T r 
is taken as October 1967, and the reference pitch angle is taken 
as 90°. It is also assumed that equation 4 can be generalized to give 

j = a r fi(L,T) E exp [-E/X r f 2 (L,T) ] (5) 


where the pitch angle dependences of the parameters are given by the 
expressions 

. m, . n 

sin (ot 0 -a c ) 

a(ao,L,T) = a r fj * a r for <f> > ao > a c 

sin (4>-a c ) 

( 6 ) 

= a r for 90° > cxq > <j> 

and 

X(ao,L,T) = X r f 2 = X r sin 11 ao/sin 11 <p for <p > «o - 

1 (7) 


= X 


r 


for 90° > otg > 
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In equations 6 and 7, m and n are L-dependent pitch angle parameters 
and tj) is an L-dependent limit for the pitch angle variation. The 
five model parameters ar, X r , m, n, and <f> are given in Table 1 for 
.05 intervals in L for L < 2 and at .1 L intervals at higher L values. 
A sixth variable (Jj) r is also given in this table, where 


00 

(Jl)r = J a r E exp(-E/X r ) dE = a r Xr (8) 

0 

This parameter is included because the variable a r has little physical 
significance, whereas (Jl) r approximates to the total unidirectional 
electron flux at the reference condition for energies above thermal 
energy. This parameter will be referred to as the electron content 
parameter; however, it is not a measure of electron density or total 
number of electrons. The hardness parameter X r is a measure of how 
fast the spectrum falls with increasing energy; that is, the larger 
X r , the harder the spectrum. The variation of (Ji) r and X r with L is 
also shown in Figures 1 and 2 , and it can be seen that the peak hard- 
ness parameter occurs at L = 1.8 and the peak electron content occurs 
at L = 1.9. In addition. Table 1 and Figure 3 show that the pitch 
angle dependences become progressively more steep with decreasing L, 
particularly at low L values where large changes in the parameter m 
can be seen. Note that no particular significance can be attributed 
to the larger changes of m as opposed to n at low L values. In general, 
both parameters are required for adequate modeling of the data, but the 
marked steepening of the pitch angle distribution at low L values could 
have been modeled by a large change in the parameter n rather than m, 
as shown in Figure 3. The parameter <(> in Table 1 appears irregular in 
its variation with L. However, the parameter <fi/a c is monotonic with L 
and indicates that the flux changes with pitch angle occur over a de- 
creasing pitch angle range as the lower L values are approached. 
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The atmospheric cutoff field values B c assumed for the present 
model are shown in Table 1 and Figure 4. For L > 1.7, B c is assumed 
to be linear with L and is given by the relation 

B c = .16 + .06L (9) 

The atmospheric cutoff based on the Jensen and Cain 1960 magnetic 
field model with a minimum altitude of 100 km is shown in Figure 4; 
equation 9 is a reasonable approximation to this cutoff. As will be 
seen from later sections, for L > 1.7 no data were available in the 
region of the atmospheric cutoff, and the influence of equation 9 on 
the fit between the model and the data is minimal. For L < 1.7, 
however, good data were available in the region of the atmospheric 
cutoff, and the B c values given in Table 1 are essentially determined 
by the data. Figure 4 shows that the B c values for L < 1.7 are below 
the values given by the linear approximation equation 9. 

The model parameters and equations presented in this section and 
the temporal parameters discussed in later sections are included in the 
APL (A Programming Language) routine AMODEL on the Goddard Space Flight 
Center IBM 360/95. A sample of the output from the program is shown in 
Table 2. The model equatorial pitch angle distribution is shown for 
three energy intervals and for one energy threshold at an epoch of 
October 1967. The energy intervals, pitch angle range, and epoch are 
input data. 

In addition to the APL programs mentioned above, a FORTRAN program 
UNI, which is capable of generating arbitrary flux from the analytical 
model, is available. Further discussion of these programs is given 
in Section 7, and more detailed information can be obtained from the 
publication by Teague (1972). While the analytical model fits the 
available data well (Section 3), some care must be taken in deter- 
mining fluxes at energy thresholds or in energy ranges different from 
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those used to evaluate the model parameters. In particular, for a 
given energy threshold or energy range, significant differences in 
the fluxes predicted in the cutoff region result from different as- 
sumptions concerning the energy range used for the pitch angle de- 
pendence. It was determined that the best results were obtained by 
using the energy limits for the pitch angle distribution corresponding 
to the measurement that was closest to the desired energy interval. 

For example, if the unidirectional flux distribution in the interval 
150 to 250 keV were required, program UNI would use the energy range 
133 to 292 keV corresponding to data from the 0G0 spectrometer (Sec- 
tion 3) for evaluation of the pitch angle dependence. This adjustment 
is performed automatically by program UNI. 

Three-dimensional SD-4060 plots of the analytical quiet day model 
flux as a function of B and L are given in Figures 5 through 7 for 
energy thresholds 


E T = 40, 100, and 250 keV 

for an epoch of October 1967. It should be noted that, while the ana- 
lytical model is applicable for electrons with energies E < 690 keV, 
the ability of this model to represent fluxes of particles above a 
given threshold, Ej, becomes suspect for Et > 500 keV for L < 1.6 ap- 
proximately, and for L < 1.4 the technique described in this section 
for determining Jj and X becomes more complex because of the signifi- 
cant hardening of the spectrum for Ej > 690 keV. This high-energy 
spectrum is discussed in detail in a later section. 

Similar plots are presented in Figures 8 and 9 for the parameters 
Jj and X. These diagrams are intended as a pictorial representation 
of the inner belt. The steepening of the pitch angle distribution at 
low L values can be clearly seen. It is also evident that the peak 
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integral electron flux above thermal energies in the inner belt occurs 
at L ~ 1.8 with a sharp falloff at low L values and a relatively slower 
falloff at higher L values. It is apparent that the model electron 
content parameter Ji has the same characteristics as the model flux in 
contrast to the hardness parameter X, which exhibits a greater falloff 
at higher L than at lower L. 
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3. DATA ANALYSIS 


For the work described in this document, data from five satellites 
were used: OGO 1, OGO 3, 1963-38C, 0V3-3, and Explorer 26. The prime 

data sets for the quiet day flux described in the previous section were 
obtained from the first three of these satellites. At a late date, 
however, data from the OV3-3 satellite became available and were incor- 
porated into the flux model with a minimum of data analysis . All of 
the satellite data were used to derive the residual Starfish electron 
model described in Section 6B. 


Data from the University of Minnesota electron spectrometers (Prin- 
cipal Investigator - Prof. John Winckler) carried on board the OGO 1 and 
3 satellites were used (NSSDC data sets 64-054A-21A and 66-049A-22A) . 
Measurements from these satellites extended over the period September 
1964 to December 1967 (Figure 10). A complete description of the 
spectrometers and the derivation of the calibration constants and energy 
ranges of the five channels have been given elsewhere (Teague, 1970). 

In this document comparison is made between the present OGO calibration 
constants and those originally presented by Pfitzer (1968), The rele- 
vant calibration constants are summarized in Table 3. The 1963-38C 
satellite was launched into an 1100-km circular orbit on September 28, 
1963, and provided data for the present model through 1967 (Figure 10). 
Data from the Applied Physics Laboratory integral electron spectrometer 
(Principal Investigators - Drs. C. 0. Bostrom and D. J. Williams) were 
used (Beall, 1969). Total particle unidirectional flux measurements 
were made, corresponding to nominal electron thresholds of .28, 1.2, 

2.4, and 3.6 MeV. Useful data for the present requirement were ob- 
tained only from the lowest energy channel because the remaining 
channels had fallen to the proton background for the time interval under 
consideration. However, the present Starfish model is compared with 
results based on data including those from the higher energy channels. 
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Calibration information for the two lower energy channels was obtained 
by integration of the efficiency curves using program EPSBART (Teague, 
1972). The results are shown in Figures 11 and 12. In each case the 
parameter e = 1/c x geometric factor is averaged over a range of 
spectral hardness appropriate to the observed spectra. Eg, Xj, and 
X 2 are the hardness parameters for an exponential spectrum, an energy 
times exponential spectrum, and an energy squared times exponential 
spectrum, respectively. For the higher energy channel the exponential 
spectrum is appropriate, whereas for the lower energy channel the energy 
times exponential spectrum is chosen (Teague, 1970). The present cali- 
bration values are also shown in Table 3 and are compared with the values 
given by previous authors. Data over the time period mid-1966 to late 
1967 were used in the present study. Explorer 26 data from the Univer- 
sity of California at San Diego detectors (Principal Investigator - Prof. 
C. E. Mcllwain) were used for the time interval January to June 1965. 

Data from two integral detectors with thresholds at 500 keV and S MeV 
were available. These data proved to be of very limited use, however, 
because it was clear that the data were in disagreement with other data 
sets. Since the data appeared less inconsistent with other data sets 
in the outer electron zone study (Singley and Vette, 1972), it is likely 
that the proton contribution to the total flux was not accurately de- 
termined. The 5-MeV data were not used in the present analysis, and 
the 500 -keV data were used only to estimate Starfish decay times over 
a limited L range (Section 6B) . EPSBART runs were also made for this 
detector, and the results are shown in Figure 13 and Table 3. 0V3-3 

data from the Aerospace Corporation differential nine-channel electron 
spectrometer were used. These data were kindly supplied by the Princi- 
pal Investigator, Dr. A. Vampola, at a late date in the evolution of the 
model and were incorporated in pnly an elementary way. The spectrometer 
provided unidirectional differential total particle fluxes in nine 
energy intervals corresponding (approximately) to electron energies in 
the range 2310 > E / (keV) > 300 (Table 3) . Data from only the two lowest 
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energy intervals were used in determining the analytical quiet day model, 
corresponding to electron energies in the ranges 375 > E/(keV) > 225 and 
600 > E/(keV) > 350. Data from the seven higher channels were used as 
the sole basis for the high-energy model discussed in Section 4. No 
separate evaluation was made of the detector calibration constants. 

As noted previously, the data sets used for the development of 
the analytical inner belt model were 0G0 1, OGO 3, 1963 -38C, and 
OV3-3. Table 3 indicates that the OGO spectrometer provided unidi- 
rectional electron flux measurements in five energy ranges . The OGO 
data sets were averaged into monthly intervals over the period 
September 1964 to December 1967 and into 0.1-L intervals over the 
range 2.0 > L > 1.3 and 0,2-L intervals over the range 2.4 > L > 2.0. 

No OGO data were available over the period October 1965 to June 1966. 

By using the first adiabatic invariant, the data were transformed into 
the form of an equatorial pitch angle distribution. The pitch angle 
coverage obtained for the time interval and L value range of interest 
is summarized in Table 4. It is apparent that the coverage obtained 
at high L values at the model epoch is poor and also that no data were 
available at this time for L < 1.5. At earlier times, however, pitch 
angle coverage was considerably better, and, with the assumption that 
the shape of the equatorial pitch angle distribution is time indepen- 
dent, these earlier data can be reduced to a common epoch for the de- 
velopment of the model. However, it is first necessary to establish 
which periods in the data interval are significantly influenced by re- 
sidual Starfish electrons and storm time electrons and which periods 
are free of these effects and are therefore quiet periods. These quiet 
periods were primarily established by investigation of the OGO data sets, 
but confirmation of these periods was obtained by analyzing the 1963-38C 
data and the 0V3-3 data. This analysis was previously described by 
Teague and Vette (1971), and a summary table of the quiet periods 
obtained is presented in Table 5. The monthly averaged data sets from 
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these quiet periods were normalized to an epoch of October 1967 using 
an cxq- independent factor. The average data sets so obtained formed 
the major basis for the analytical model. The monthly averaged 0G0 data 
sets indicated that the counts for energies greater than 690 keV regis- 
tered nonzero values only for periods when Starfish and storm time 
effects were significant. It was therefore decided to develop a quiet 
day model valid for energies E < 690 keV. Observations from the 0V3-3 
data indicated significant fluxes for E > 690 keV, however, and it 
became apparent that the OGO satellite generally was not measuring 
these because of a combination of a high intensity threshold and the 
sampling of only low equatorial pitch angles. 

The 1963-38C data were available in the form of total counts as a 
function of time for a specified narrow B ± interval C - 02 to .01 gauss) 
at .05 intervals in L for 1.2 < L < 1.6 and .1 intervals for 1.6 < L 
<2.0. In addition, perpendicular flux distributions were available 
at certain epochs over the period 1966 to 1967. The approximate equa- 
torial pitch angle range covered by these distributions is indicated in 
Table 4. The proton background was removed from the total flux measure- 
ments using the 1963-38C measurements given by Beall (1969). In general, 
flux distributions were not available for quiet periods, and it was 
assumed that the shape of the storm time distributions was identical to 
that on a quiet day. Average perpendicular flux distributions were 
obtained by normalizing to epoch October 1967 using count rates as a 
function of time for a fixed B^ interval. 

As noted previously, the 0V3-3 data did not become available until 
much of the modeling activity had been performed. Dr. Vampola pro- 
vided microfilm plots of equatorial perpendicular flux as a function 
of time for the period late 1966 through 1967 for each of the nine 
channels. Complete data analysis had been performed by Dr. Vampola. 

The unidirectional flux measurements after removal of the proton 
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background had been reduced to equatorial perpendicular fluxes by use 
of pitch angle distributions based on the 0V3-3 data. For L > 2, how- 
ever, Dr. Vampola had made the assumption of linear pitch angle depen- 
dence owing to poor data coverage. The assumed pitch angle dependences 
were not available to NSSDC at the time of the modeling activity. Only 
equatorial perpendicular flux values for the two lower energy channels 
(Table 3) were used for the analytical quiet day model . 

Plots of the 0V3-3 data used in developing the present quiet day 
model are given in Figures 14 and 15, and CalComp plots of the 0G0 and 
1963-38C data are presented in Figures 16 through 57. The averaged 
count rate distributions given by the 1963-38C data are shown with the 
equatorial pitch angle distributions from the 0G0 data for the quiet 
periods previously discussed. The OGO data sets are not shown normalized 
to a common epoch in order to illustrate temporal flux changes. These 
changes are the subject of following sections. The ordinate scale on 
Figures 16 through 57 is arbitrary count rate. Pfitzer (1968) performed 
normalization of the absolute count rates measured by OGO 1 and OGO 3 in 
order to obtain a common data base. Furthermore, the 1963-38C data 
shown in these figures do not provide absolute count rates because they 
are normalized to a common epoch. 

The 0V3-3 data are shown as a function of L in Figures 14 and 15 
for the mid-1966 and early 1967 periods, respectively. To reduce the 
model to the correct epoch, the flux ratios determined from the OGO 3 
690 > E/(keV) > 292 data were used. In general, more data were avail- 
able for the earlier epoch shown in Figure 14, and more emphasis is 
given to these data in determining the model parameters. The broken 
line in Figure 14 represents the estimated background in August 1966, 
and the full line indicates this background with the Starfish residual 
flux added (Section 6B). In Figure 15 the broken line represents the 
model for epoch August 1966, and the full line represents the model 
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for epoch April 1967. A small Starfish residual not shown remains at 
low L values for epoch April 1967. However, the data at these L values 
remain slightly above the model. It can be seen that the agreement be- 
tween model and data in Figure 15 is poorer than in Figure 14, with the 
data usually below the model. As noted, the model parameters were de- 
termined giving greater emphasis to the August 1966 epoch because of the 
more limited data available at the April 1967 epoch. However, the com- 
parison shown in Figure 15 does cast some doubt on the accuracy of the 
solar cycle parameters obtained with the 0G0 data for high L values at 
these energies (Sections 6A and 7) . 

In general, the data from the OGO 3, 1963-38C, and 0V3-3 satellites 
provide reasonable pitch angle coverage, as can be seen from the tables. 
However, in a number of regions data are poor, and the analytical model 
is no better than extrapolation from higher pitch angles or other L 
values. For example, at high L value (L > 2 approximately), poor pitch 
angle coverage at the model epoch is obtained from the OGO data, as 
shown in Figures 48 through 57. This deficit is in part offset by the 
better coverage at earlier quiet periods and by the equatorial perpen- 
dicular measurements from the 0V3-3 data. That is, the shape of the 
pitch angle dependence is well deteimined, but the absolute fluxes at 
the model epoch may be suspect. Further, at low L values (L < 1.4 ap- 
proximately) the OGO data are no longer useful, and thus limited spectral 
information is available. Fortunately, at these L values excellent pitch 
angle coverage is afforded by the 1963-38C data, although the background 
proton flux becomes significant in comparison to the electron flux, and 
the data standard deviation increases accordingly. Again there is the 
problem of reducing the 0V3-3 and the 1963-38C data to a common epoch, 
which results in some uncertainty in the absolute value of the electron 
content parameter (Ji) r . At intermediate L values where the coverage of 
the independent variables is best, some conflict arises between the 
1963-38C data and the OGO data in the region of the atmospheric cutoff, 
thus resulting in some uncertainty in the pitch angle parameters m and n. 
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A case in point is L = 1.6, Figures 32 through 35. It is apparent 
from the comparison of the 0G0 data and the model plots that the model 
plots indicate a more rapid fall at low equatorial pitch angles than 
the 0G0 data would indicate. In Figure 35, however, the discrepancy 
between the 1963-38C data and the model is seen to be the reverse; 
i.e., the model is slightly higher than the data would indicate in 
the region of the atmospheric cutoff. The conflict between the two 
data sets occurs at L values in the range 1.9 > L 2 1.6, and the 
1963-38C data are generally considered to be more reliable. The model 
is primarily a fit to these data for low L values. It was noted in 
the previous section that at low L values the atmospheric cutoff is 
essentially determined by the 1963-38C data and that B c in this region 
deviates from the linear assumption for higher L values . Very signifi- 
cant model flux changes occur in the cutoff region as small changes in 
B c are made, owing to the (ao~a c ) term in equation 6. Since this region 
is ill defined at high L values, it is probable that the model flux is 
associated with a considerable error for low pitch angles at these L 
values. An error of a factor of 3 or 4 is not unlikely. 

The equation for the differential spectrum, equation 4, was chosen 
on the basis of the work performed by Teague and Vette (1971). A 
number of functional forms describing the pitch angle dependence of 
the spectral parameters were tested. Equations 6 and 7 were chosen 
because they best describe the flux change with pitch angle. Other 
functions can be used that better describe the pitch angle dependence 
at low or high L, but these equations represent the best compromise 
for the total range. In addition, the inclusion of the a c parameter 
in equation 6 ensured that the model flux became zero at oq = a c . 

The model parameters were fitted to the data by hand iteration using 
the APL program AMODEL (Teague, 1972) previously discussed. Also, a 
separate APL program, FLUXOV (Teague, 1972), giving the equatorial 
perpendicular flux was written for the OV3-3 data. An integrated form 
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of the differential spectrum, equation 4, was used for fitting to the 
measured fluxes in a given energy interval . That is , the average dif- 
ferential flux, ji, in the energy interval E ±2 to Eii is 


(J T ) 

I-'r 

E. ,,-E. . 
i2 ll 



Fits were initially performed in the areas where the data coverage was 
good, i.e., intermediate L values. The results for these L values were 
used to establish trends for estimating the model parameters in L 
regions where the data coverage was poor. At low L values for which no 
0G0 data were available, the least squares APL library program NONLIN 
(Wilson, 1971) was used for fits to the 1963-38C data. The 0V3-3 data 
were used exclusively for determining the cutoff $ for the pitch angle 
functions. In addition, the data set provided valuable confirmation 
of the model parameters determined from the OGO data. A complete 
description of the APL programs mentioned here can be found in the 
document "Inner Zone Computer Programs" (Teague, 1972). 

Figures 14 through 57 show that in general the agreement between 
the model and the data is good, and it is considered that the model 
flux is within the standard deviation of the data. It is further 
considered that the model flux provides as good an estimate of the 
data at high energy (E a 690 keV) as at low energy; that is, the 
spectral function is not a simply linearly weighted fit to the data. 

In Figure 15 agreement between the model and the low-energy 0V3-3 
data is not particularly good. However, the data have a large standard 
deviation, and the model is reduced to an epoch of August 1966 by using 
B-independent factors based on OGO 3 data that are associated with a 
significant error. This point is discussed further in a following 
section. 


18 



Pfitzer (1968) performed polynomial least squares fits to the 0G0 
data sets to model the pitch angle dependence. It is interesting to 
compare these polynomial fits with the present model. Pfitzer's fits 
were performed to both 0G0 1 and OGO 3 data sets, and a B-independent 
factor was used to normalize the different epochs , The comparison with 
the AE-5 model for October 1967 is made by normalizing at the equator, 
and the results are shown in Figures 58 and 59 . It should be noted that 
differences between the two pitch angle distributions at low ccq may be 
artificial owing to the normalizing process. The slopes can be compared, 
however, and it is apparent that the present pitch angle dependence gives 
a greater slope for most cases. This result is partly due to the inclu- 
sion of the 1963-38C data in the present model, as noted previously. 
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4. HIGH-ENERGY ELECTRONS 


The preceding sections have discussed the development of a quiet 
day analytical model appropriate for electrons with E < 690 keV. As 
noted in Section 3, a limited amount of data from OV3-3 were available 
for electrons in the energy range 2.31 > E/ (MeV) > 0.3. It is apparent 
from Table 3 that the OGO spectrometer sampled the flux for E > 690 keV 
in the two higher channels. However, for those periods not influenced 
by Starfish electrons (Section 6B) , these channels generally registered 
no flux because of a combination of high intensity threshold and low 
equatorial pitch angle sampling. In this section, the development of 
the OV3-3 data into a high-energy model is discussed. It should be 
emphasized that this portion of AE-5 is based on only a single data set 
that provides no information concerning the pitch angle distribution. 
This portion is therefore somewhat provisional, and development of a 
more comprehensive model for high energies has already begun using ad- 
ditional data from other experiments that have become available. 

In addition to the lack of data, a fundamental problem of electron 
modeling arises for the high-energy electrons in the inner belt. The 
inclusion of the effects of magnetic storms in an electron model is 
necessary for those regions of B-L-E space where these effects signif- 
icantly influence the average environment. Three variables must be 
considered in determining the influence of magnetic storms : (1) the 

frequency, (2) the flux change in relation to the undisturbed (quiet 
day) background, and (3) the duration. Assessment of the importance of 
magnetic storms can be performed in practice with consideration of the 
first two variables alone because these exhibit much greater variation 
with E and L than does the third variable. 
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The implicit assumption of the quiet day analytical model presented 
in Section 2 is that storm effects do not significantly affect the aver- 
age environment for E < 690 keV and 2.4 > L > 1.2. It will be shown in 
this section that the frequency and magnitude of the flux change caused 
by magnetic storms are such that this is a true statement with the ex- 
ception of the higher L values at certain energies. As L increases 
through the slot region and into the outer belt, both the frequency and 
relative flux change due to storms increase rapidly, until at L ~ 3.0 
it is impossible to isolate quiet day periods from disturbed periods. 

Thus the only modeling approach available in the outer belt is to 
average the electron fluxes, including the storm effects, and to treat 
the significant excursions of the fluxes from this average statistically. 
This is a valid technique, provided that the averaging period is suf- 
ficently large to represent a random sample of the storm effects, because 
a model user will be applying the environment to a period other than the 
averaging period. In general, the frequency of effective magnetic storms 
is sufficiently high in the outer belt for a random sample to be obtained 
with a relatively short averaging period. The latest outer belt model 
AE-4 averaged over the period mid-1966 to late 1967 for this purpose. 

In principle, variations of storm effects with solar cycle can be in- 
cluded in both the average value and in the standard deviation, although 
in practice for AE-4 it has been possible to derive only average flux 
conditions approximating to solar maximum and solar minimum owing to 
the data coverage (Singley and Vette, 1972). 

As opposed to the outer belt, the frequency of effective magnetic 
storms in the inner belt is too low for a statistical approach. Over 
the period June 1966 to December 1967, for example, excursions of the 
flux above the quiet day level were observed for approximately 20% 
of the period. The problem of what approach to adopt arises when the 
frequency of effective magnetic storms is low and the relative flux 
change is high. It will be shown that this situation arises for high- 
energy electrons in the inner belt. In these circumstances the flux 
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varies considerably from quiet to storm condition and in such a way that 
the changes from one condition to the other are unpredictable and cannot 
easily be modeled. An average storm model requires an impractically 
long averaging period to obtain an effectively random sample. In addi- 
tion, before the average environment becomes meaningful to the user, he 
must fly his satellite in that environment for a similarly long period. 
The most important variable becomes whether the satellite is flying in 
the quiet or the storm environment. The most meaningful approach may 
be to provide the user with (1) a quiet day environment and (2) the 
ability to obtain a storm environment corresponding to a storm of an 
arbitrary magnitude. That is, a storm model correlated with some mag- 
netic index or parameter would be desirable. This approach is presently 
being investigated, but, because of the lack of data available for the 
present model, the effects of storms are included in a crude averaging 
fashion, as discussed below. However, it should be appreciated that 
the present high-energy model at best provides only an approximation 
to the environment. 

The energy intervals covered by the 0V3-3 spectrometer are shown 
in Table 3. As noted previously, no calibration information has yet 
been provided, but the data were supplied in the form of equatorial 
perpendicular flux. To assess the importance of magnetic storms, quiet 
day flux levels were determined at a variety of epochs together with 
the average storm time flux based on the period June 1966 to December 
1967. It was considered that the ratio of these two fluxes was in- 
dicative of the importance of magnetic storms . 

The quiet day high-energy equatorial flux is shown in Figures 60 
and 61 for two epochs August 1966 and October 1967 for selected L 
values. Some data from April 1967 are included in epoch October 1967. 

The data are plotted at the midpoint of each energy range. Two dis- 
tinct temporal variations can be observed at different L values. In 
Figure 60 it is clear that the August 1966 data are lower than the 
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1967 data, as was observed at lower energies in the 0G0 and 1963-38C 
data. Above 1.49 MeV this temporal variation is not observed. The 
standard deviation of the data in this region becomes large, however. 

The August 1966 data for L = 2.0 and 2.2 are not shown but exhibit a 
similar temporal variation. In Figure 61, the data at L = 1.8 show a 
similar increase with time but, at L < 1.6 approximately, a reduction 
in flux with increasing time is observed. This apparent decay is in- 
terpreted as depletion of residual Starfish electrons. These temporal 
variations are discussed further in following sections. 

Integral spectra were derived from the differential fluxes, and 
the quiet day integral fluxes are shown for an epoch of October 1967 
on Figures 62 and 63 for L > 1.7 and L < 1.7, respectively. The data 
coverage at the epoch of October 1967 is poor, and a large portion of 
the data shown on Figures 62 and 63 are based on the early quiet epoch 
of April 1967. Where data are available at both epochs, any temporal 
variations are obscured by standard deviation in the data. Estimates 
were made of the flux above the energy range measured by the detectors. 
At low L values for which the spectra were hard, the contribution of 
this part of the spectrum to the integral flux above 2.31 MeV was 
significant in relation to the error of the data, which is large at 
the high-energy end of the spectrum. For L > 1.8 approximately, fluxes 
at energies above those measured by 0V3-3 did not contribute signifi- 
cantly to the integral spectra up to thresholds of 2.31 MeV. Figures 
62 and 63 show that the spectrum becomes harder at low L values and 
that the flux exhibits a maximum at L = 1.3-1. 4. That is, the radial 
profile at the high energies is markedly different from that at the 
lower energies discussed in previous sections that exhibit maxima at 
L = 1.8-1. 9. 
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To assess the importance of magnetic storms, average flux values 
were derived for the period June 1966 to December 1967, including three 
storm periods corresponding to the magnetic events beginning in late 
August 1966, early January 1967, and late May 1967. The average fluxes 
F were derived using the simplified model given in Figure 64. It was 
not possible in each case to measure all the variables shown in this 
figure, and many estimates were made from the available data, parties 
larly concerning the storm decay times. From these results the ratio 


Average Flux June 1966 to Dec. 1967 
Quiet Day Flux Oct. 1967 


(H) 


was derived. The ratio Rg is shown as a function of energy for a variety 
of L values in Figure 65 and is based on differential fluxes. This ratio 
is plotted at the midpoint of each 0V3-3 channel. It is clear conceptu- 
ally that Rg is only a qualitative index, and considerable error may be 
associated with actual R s values presented in Figure 65. It can be seen 
that the peak storm effect is observed in the energy range 1100 > E/keV 
> 850 for 2.8 > L > 1.8, with sharp falloff with both increasing and 
decreasing E. Furthermore, the storm effect increases markedly as L 
increases. Below L = 1.8, Rg was essentially unity. Of course, that 
does not imply that storm effects cannot be observed at these L values, 
but rather that their significance in the context of providing an aver- 
age model is small and within the standard deviation of the quiet day 
model . 

From Figure 65 it may be concluded that the quiet day model pre- 
sented in Section 2 for 690 > E/keV > 36 provides a good estimate of 
the environment for L < 1.9 and that the quiet day model of electrons 
with E > 690 keV presented in this section provides good flux estimates 
for L < 1.8. For L values above these ranges, the occurrence of a 
magnetic storm is likely to increase the average flux by a significant 
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amount. It is interesting to consider the probability that the observed 
flux in the inner belt will be above the quiet day level or above the 
average flux. The probability distribution is determined largely by 
the idealized storm model used (Figure 64), and the probabilities are 
heavily dependent on the storm decay times . Figure 66 shows the prob- 
ability that the flux is above the quiet day background (epoch October 
1967) and above the average flux F as a function of energy for L = 1.9 
and L = 2.8 for the period June 1966 to December 1967. Clearly, the 
probability that the flux will be above the quiet day background P(h^) 
is quite significant. Tbe probability that the flux will be above the 
average, P(F), is significantly lower, having a typical value for the 
L range 2.8 > L > 1.8 of approximately 14% for most energies. Note 
that no clear differences become apparent in P for different L values . 
Significant increases in the probability occur at low energies (L = 2.8, 
Figure 66) arising from increases in the storm decay times, but gener- 
ally these are relatively unimportant in the inner belt because they 
coincide with small values of Rg (Figure 65) . 

As a crude estimate of storm effects, the average storm time flux 
F is presented as the high-energy model. As noted above, if the period 
June 1966 to December 1967 can be regarded as typical, the probability 
of the flux being above this level is approximately 14%. However, it 
should be appreciated that this probability reduces only slowly with 
increasing flux level owing to the spiked nature of the storms. 

Figure 67 shows the averaged equatorial perpendicular integral spectrum 
for L > 1.7 and can be compared with the quiet day fluxes presented in 
Figure 62. At L = 1.7 the spectrum remains unaltered, but at higher 
L values the inclusion of storm effects results in a flattening of 
the radial profile to the extent that little change in flux occurs 
between L = 2.2 and L = 2.4, with the flux at L = 2.4 becoming 
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greater than that at L = 2.2 for 1.8 > E/(MeV) >0.7. For L < 1.7 
the model spectra are unaffected by magnetic storm effects, and the 
spectra presented in Figure 63 are compatible with those shown in 
Figure 67. 
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5. THE OMNIDIRECTIONAL FLUX MODEL AE-5 


In this section the development of the complete model in the form 
of omnidirectional flux and based on the information presented in 
Sections 2 and 4 is discussed. This model, designated AE-S, is valid 
for the epoch October 1967, for 1.2 < L < 2.8, and for 4.0 > E/(MeV) 

> .04. Efforts are made to ensure that the inner zone model AE-5 is 
compatible with the latest outer zone model AE-4 (Singley and Vette, 
1972) at the interface slot region L ~ 2.6. 

A. The Inner Zone Model 

The basic constituents of this model have been discussed in Section 
2 Cthe analytical model for E < 690 keV) and in Section 4 (the high- 
energy electron model. The form of the model is omnidirectional in- 
tegral flux provided as a function of B, L and energy threshold E T . 

The effects of magnetic storms and the high-energy electrons were dis- 
cussed previously. It is clear from Figure 65, however, that storm 
effects are important at lower energies for L > 1.9, and thus some 
modification to the quiet day analytical model is necessary. In addi- 
tion, the ability of the analytical model to estimate integral flux 
becomes suspect because the portion of the spectrum above 690 keV is 
estimated by fits to data for E < 690 keV. Figure 68 shows the quiet 
day and the average storm time models for L > 1.9 and 900 > E/(keV) > 

200. In this figure the analytical model is used to give the quiet 
day integral flux for E-p < 500 keV, approximately. Above this energy 
the quiet day spectra are matched with those given by the high-energy 
0V3-3 data. In general, only minor modifications have been made, that 
is, within the standard deviation of the analytical model and the 0V3-3 
data. The full lines in Figure 68 represent the final equatorial 
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perpendicular integral flux model, and the modifications to the integral 
flux for E < 500 keV to include the average storm effects can be assessed 
by comparison with the quiet day flux lines. In general, the analytic 
model can be used for determining integral fluxes up to energies Ej in 
the range 300 to 400 keV dependent on the L value. 

At L < 1.9 storm effects become insignificant. However, the quiet 
day analytic model requires significant modification because of the 
high-energy portion of the spectrum, which becomes increasingly hard 
at low L values. The integral flux spectra are indicated in Figure 
69 for L = 1.8 and 1.7. In each case the analytic model is used for 
integral flux up to Ej - 500 keV. It can be seen that no modification 
is made at L » 1.7 for either storm effects or high-energy electrons, 
and only minor storm effects are evident at L = 1.8. It is interesting 
to note the crossover of the spectra at 560 keV that also occurs at the 
lower L values shown in Figure 70. In this figure it can be seen that 
considerable correction is required to the integral flux for E T < 500 
keV owing to the hardening of the spectra at high energy. For example, 
at L = 1.3 the integral flux at 500 keV is increased by an order of 
magnitude because of this effect. Despite these large corrections, 
there is no conflict between the differential spectrum given by the 
analytic quiet day model at 500 keV and the 0V3-3 data. This is il- 
lustrated in Figure 71 for L = 1.3 through 1.5. In Figure 71 it is 
apparent that the effects of the high-energy electrons for L < 1.4 are 
such that the method used for determining the quiet day model (Sections 
2 and 3) is not applicable for these L values. This is because the 
average differential flux measurements of the 0G0 satellite are not 
compatible with the integral flux measurements of the 1963-38C satel- 
lite. Further, for L < 1.3, estimates of the quiet day parameters Jj 
and X cannot be obtained from the 1963-38C data. To account for the 
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effects of high-energy electrons, an iterative procedure is adopted, 
and Figures 20, 21, and 25 indicate the predictions of the quiet day 
model based on Jj and X and the final model, including the effects of 
the high-energy electrons . 

To make the present AE-5 model compatible with the outer zone 
model AE-4 (Singley and Vette, 1972), the high-energy data are extra- 
polated to provide flux estimates for particles with energies up to 
= 5 MeV. The extrapolations performed are shown in Figures 72 
and 73 for L > 1.6 and L < 1.6, respectively, The extrapolation is 
performed until the flux has fallen to 10 -1 electrons/cm 2 -sec-ster 
(Section 6). The outer zone electron model AE-4 provides insignifi- 
cant fluxes at all L values for E T > 4,85 approximately. The present 
extrapolation is performed with a similar criterion. Evidently the 
extrapolation procedure performed above leads to a considerable error 
for E T > 3 MeV, and the error at 4 MeV may be as much as an order of 
magnitude. 

The analytic quiet day model provides the pitch angle dependence 
for electrons with E < 690 keV. It is assumed that the storm contri- 
bution added to this model at high L values does not change this pitch 
angle distribution. Examination of the storm time pitch angle distri- 
bution from the 0G0 data supports this assumption. As noted in Section 
3, no pitch angle information is available for E > 690 keV. It is 
assumed that the observed pitch angle dependence for the energy range 
690 > E/(keV) > 292 can be extended to higher energies. This assump- 
tion is supported by the fact that the observed pitch angle dependence 
in the outer zone is energy independent. This assumption is likely to 
result in some error, however, since the pitch angle distribution in 
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the inner radiation zone shows variation with energy for E < 690 keV. 
Three-dimensional SD-4060 plots of the high-energy model flux as a 
function of B and L are given in Figures 74 through 77 for energy 
thresholds 


Ej - 0.5, 1.0, 2,0, and 4 . 0 MeV 

Using the equatorial perpendicular flux values presented in this 
section and the analytic model pitch angle distributions, omnidirec- 
tional integral fluxes J were obtained performing the integration 

B C 

JCB.L.Et) = j (12) 

B 

where is the perpendicular flux corresponding to B x . The resulting 
inner belt radial profiles are presented in Figures 78 and 79 for 
selected energies. It is clear that the peak flux moves to lower L 
values as the energy increases. The hardening of the high-energy spec- 
trum at low L value results in a doubly peaked profile at Ex = 500 keV. 
This profile would not be expected from simple diffusion mechanisms, 
and it is considered that there is a significant Starfish contribution 
to the flux at these energies at low L values for an epoch of October 
1967. This subject is discussed further in Section 6B. 


B. The Interface Region L ~ 2.6 

The above inner belt model is valid for 1.2 < L £ 2.4. Model 
AE-5 is extended to include L = 2.6 earth radii by comparison with the 
outer zone model AE-4 (Singley and Vette, 1972). No data are considered 
at this L value, and the model is determined solely by interpolation 
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between L = 2.4 and L = 2.8 or 3.0. The radial profile in the inter- 
face region is shown in Figures 80 and 81 for E < 1 MeV and E > 2 MeV, 
respectively. Figure 80 shows that the two models provide similar 
flux levels on either side of the interface region and that the trends 
given by the radial profiles are compatible. The minimum flux at any 
energy is observed to move to lower L values with increasing energy. 
Figure 81 shows that there are considerable differences between the 
flux levels on either side of the interface region at the higher 
energies. That is, the outer belt spectrum is significantly harder 
than the inner belt spectrum. However, we note that there is con- 
siderable error associated with the inner belt spectrum above 3 MeV 
because of the extrapolation process discussed in Section 5A. The 
trends given by the radial profiles of the two models are apparently 
contradictory for E T > 3 MeV. The outer zone model at L = 2.8 at 
these energies is based entirely on Explorer 26 data and, as noted 
in Section 3, these data are anomalous at L values where the proton 
background becomes significant. In the final model shown in Figure 
81, the outer zone model is modified at L = 2.8 for Ej > 3 MeV to 
make it compatible with the inner zone model. The interpolated 
spectrum at L = 2.6 is shown in Figure 82 and the B dependence is 
shown in Figure 83. In Figure 83 the fluxes are normalized to 
unity at the equator. 


33 



6. TEMPORAL VARIATIONS IN THE INNER RADIATION BELT 


Isolation of the periods remote from the influence of magnetic 
storms has been discussed in previous sections for energies E < 700 keV 
approximately. Data from these periods exhibit temporal behavior at- 
tributable to decay of Starfish electrons and solar cycle effect. Sep- 
aration of these two effects is not a simple process because a large 
portion of the data is significantly influenced by both, which results 
in an iterative procedure that also affects the model parameters de- 
scribed in the previous sections. In some cases, however, data that 
are predominantly affected by one or other of the variations can be 
isolated as in late 1964, when Starfish decay dominated for L < 2.0, 
and as in mid- 1966 and later, when solar cycle variation dominated 
for low energies . The two temporal variations are discussed separately 
in the following sections with considerable cross-referencing between 
sections . 

A. Variations with Solar Cycle 

The analysis in this section is concerned primarily with electrons 
of E < 690 keV. It will be seen that similar variations can be detected 
at higher energies, but in many cases in the context of an electron 
model these variations are small in comparison to the effects of mag- 
netic storms. 

It is assumed that the solar cycle effect is not B dependent. Then, 
the ratio of the unidirectional flux at the epoch of October 1967 to 
that at time T months measured from solar minimum (assumed to be 
September 1964) is defined as R<p (E,L). It is apparent in some instances, 
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however, that the solar cycle effect is B dependent, with the flux at 
the lower pitch angles exhibiting less variation than at higher pitch 
angles (e.g., Figures 42 and 46). In general, however, insufficient 
data are available for modeling this variation, and B independence is 
assumed. 

It can be seen that at L = 2.4 (Figures 55 through 57) Starfish 
electrons with E s 690 keV have decayed to insignificant levels by 
early 1965, and consequently both the 0G0 1 and 0G0 3 data sets can be 
used for investigation of the solar cycle variations. Figure 84 shows 
the model ratio Rj for this L value. From Figures 55 through 57 it can 
be seen that the data are poor for L = 2.4 and T = 37 (October 1967), 
and the model for this time is considerably influenced by lower L values 
and by earlier data for which the pitch angle coverage is markedly 
better. It should be emphasized that there is very significant error 
associated with determining the ratios Rj that is caused primarily by 
standard deviation in the data and secondarily by unmodeled B dependence, 
and thus care must be taken in attributing significance to the shape of 
the curves presented in Figure 84 and to those subsequently presented 
for the variable Rj. In extreme cases this error may be 50% for energies 
690 > E/(keV) > 292 and larger for lower energies where data are poor 
and Rf is smaller. From Figure 84, however, it can be concluded that 
the solar cycle effect increases markedly with energy up to 690 keV, 
and later it is shown that the L-dependent characteristics can be ob- 
served. In addition, in Figure 84 the broken lines indicate the ratio 
R? for linear variations of flux with time at L = 2.4 based on the first 
data point in each case at T = 27. At earlier times it is clear that 
for the two highest energy ranges the observed ratios Ry are signifi- 
cantly lower than the linear values, thus indicating a reduction of 
the solar cycle effect towards solar minimum, as would be expected. 

The reverse trend is observed for the lowest energy, and although it 
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is repeated at lower L values , it may simply result from error in the 
ratio Rt at T = 22. However, it is evident that the solar cycle effect 
is approximately linear for 22 < T < 37 independent of energy. 

Similar qualitative observations can be made for L = 2.0 and 2.2 
because Starfish electrons cease to be significant by mid-1966. Further, 
at these L values the pitch angle coverage is more comprehensive and more 
accurate values of Rp can be determined. For L < 2.0, however, the flux 
at T = 22 becomes influenced by Starfish electrons in the energy range 
690 > E/(keV) > 292. For the energy interval 292 > E/ (keV) > 133, Star- 
fish electrons cease to be significant by mid-1965 (Section 6B) and all 
the OGO 3 data can be used for assessing the solar cycle variation. For 
the lowest energy range. Starfish electrons have decayed to insignifi- 
cance by late 1964, and both the OGO 1 and OGO 3 data sets can be used 
for solar cycle effect. For the high-energy data containing a signifi- 
cant Starfish contribution, estimates of this contribution can be made 
using the analysis presented in Section 6B. As noted, these estimates 
are also a factor in determining the quiet day model parameters dis- 
cussed in previous sections. The resulting quiet day flux distributions 
are shown in Figures 14 through 57, labeled "T = 22 EST BACK" together 
with the total model flux and the data at T = 22. This estimated back- 
ground can be used to determine an approximate Rp value for T = 22. In 
addition, for the highest energy channel estimates of Rp at T = 31 and 
27 can be made using quiet day data from these epochs. As noted pre- 
viously, for L = 2.4 the flux varies approximately linearly with time 
for 22 < T < 27 for all energies. The L dependence of the ratio Rp is 
shown in Figure 85 for various epochs for the three channels of OGO data 
and at T = 22 for the E > 255 keV 1963-38C data. In general, Rp in- 
creases with L and reaches a peak at L = 2.0 approximately, with a 
sharper falloff at lower L than at higher. The peak at L = 1.7 in the 
133 > E/(keV) > 36 is not considered significant owing to the standard 
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deviation in the data. A summary of the ratios Rj determined from the 
OGO 3 and 1963-38C data is shown in Table 6 for T = 22 months. In 
addition, carpet plots of Rt(T,L) are given for energy ranges corre- 
sponding to the OGO satellite in Figures 86 through 88. (Note that a 
brief description of carpet plots is given in Appendix A.) Some 
smoothing of the data presented in Figures 86 through 88 and Table 6 
has been performed, particularly at low energies and low L values. It 
is considered that the resulting changes are small in relation to the 
standard deviation of R T . In determining the carpet plots for Rt, the 
Starfish model described in the next section has been used to estimate 
the contribution from Starfish electrons to the total flux at various 
values of T. As noted above and shown in Figures 14 through 57, this 
includes L < 1.9 and T < 22 for 690 > E/(keV) > 292 electrons. In ad- 
dition, the same process is performed for 292 > E/(keV) > 133 electrons 
at T = 6. Where the Starfish contribution to the total flux is large, 
the possible error in the ratio Rj becomes very large, and the early 
epochs with low L value shown in Figures 87 and 88 must be regarded as 
provisional. However, it is interesting that a certain consistency exists 
in the values of R t at early epochs. For electrons 690 > E/ (keV) > 292 
approximately, 85% of the total increase in flux between T = 6 and 37 
occurs over the period T = 22 to 37. This approximate figure of 85% is 
independent of L. A corresponding effect is observed for lower energies 
with 75% for 292 > E/(keV) > 133 and 27% for 133 > E/(keV) > 36. 

E/(keV) > 36. 

To provide solar cycle parameters compatible with the model pre- 
sented in previous sections, it is necessary to derive values of Rj. 
as functions of a given energy threshold E^. These values can be 
based only on those presented in Figures 86 through 88 and determined 
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from the OGO data. The integral flux at an arbitrary" epoch T can be 
expressed in terms of the OGO energy ranges as 


k 

jCT,E T ) = l j-CE.^E.^T) + j (>690 ,T) 


(13) 


where E^j = Ej, Ej^ - 690 keV, and Ej^ is the upper energy limit of the 
channel measuring electrons with energy E = Ej. From Figures 60 and 61 
it is clear that particles with E > 690 and L > 1.6 exhibit a solar 
cycle effect that decreases with increasing energy and also that Star- 
fish decay is evident at lower L values. It is also clear from the 
discussion of magnetic storms, however, that this solar cycle effect 
is unimportant in relation to these storm effects, and for the present 
model it has not been possible to model the variation of these effects 
with solar cycle. The effect of the Starfish decay process at low L 
values is also unknown, and it is therefore assumed that j (>690) is not 
dependent on time. In general, however, this term has a relatively 
small effect on the value of R-j. for integral energy fluxes. In addition 
to storm particles with E > 690 keV, it has been shown that significant 
storm contributions exist at lower energies for high L values. It is 
assumed that these contributions are also time independent and that 
equation 13 can be written as 



k 


k 

11 

/ — \ 
E-< 
W 

V * 

♦t“v 

l JiCT.E^) 

+ 

q 



s + j (>690) 


(14) 


o 


39 



In general, a solar cycle ratio Rj (Ej) for a given energy threshold can 
be given by the expression 


j(37,E T ) 

j(T,E T ) = Rt( - t>Et ) 
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(15) 

where J(37 ,Et) is given by equation 13 with T = 37, and (RtH is the 
solar cycle ratio for Channel i presented in Figures 86 through 88. 

There are obvious errors involved with equation 15, particularly for 
regions of E,L space where significant modification to the quiet day 
spectra is made by high-energy or storm time particles that are assumed 
independent of time. Equation 15 can be used to relate to the 0G0 and 
1963-38C solar cycle parameters presented in Table 6. In practice 
simpler forms of equation 15 can be given for certain L value regions. 

At L values where the storm contribution is insignificant below 690 keV, 
Rt(Et) becomes 

R t CT,E t ) = j(37,E T ) 
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and at L values for which storm effects become most significant as 


R t (T,E t ) = j(37,E T ); 
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The comparison between 1963-38C and the predictions of equation 15 
is shown in Figure 89. The agreement is well within the standard devia- 
tion of 1963-38C data, which is typically 25%. The term [j CE-jOJg in 
equation 17 was not a factor in the present calculation, although 
j 0*690) in equation 16 became significant for all L values. Equations 
16 and 17 were used to derive carpet plots of RjC e t» t ) E T = 40, 1°0, 
250, and 500 keV, These plots are presented in Figures 90 to 93. Note 
that the carpet plot technique could not be employed for Ep = 500 keV 
because of the shape of Rp(T) at the various L values. For Ej > 250 keV, 
the storm term is unimportant, but the high-energy term j (>690) is sig- 
nificant, particularly at low L values. At 500 keV both terms became 
very significant, especially at low and high L values, and Figure 93 
must be regarded as very provisional because of unmodeled temporal 
variations of these terms. In general the effect of these terms becomes 
most important at low values of T, and thus the early epochs are less 
reliable than the later epochs. 

This section has been concerned with determining solar cycle ratios 
for T < 37 months. A more important and less well-defined problem, in 
the context of the usage that the present electron model will receive, 
is extrapolation for T > 37 months. Of course, there is no guarantee 
that solar cycle ratios determined for one particular solar cycle are 
applicable to another. Also, owing to asymmetries in the solar cycle 
about solar maximum, error is likely to be introduced by the assumption 
that ratios based on data between solar minimum and maximum (as with the 
present values) can be applied to a time period between the maximum and 
the subsequent minimum. However, the errors due to these sources are 
small in comparison with the considerable error associated with the 
present values of Rj. These values of R-p extend over the period 
37 > T > 22 ; that is, they do not cover the complete period solar mini- 
mum to solar maximum. Figure 94 shows the time-averaged sunspot number 
1^ as a function of time for the present solar cycle. It is evident 
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that the present epoch of T = 37 coincides approximately with the 
maximum R z observed for this cycle. Further, the earliest value of Rj 
at T * 6 coincides approximately with a minimum in R z . Therefore, 
coverage of the period 37 > T > 22 incorporates the majority of the 
flux increase due to solar cycle for the present cycle. From Figure 
94 it is clear that R z remains essentially constant over the period 
1968 to mid-1970. If the inner belt solar cycle effect is well-corre- 
lated with sunspot number, it may be assumed that the flux remains 
constant at the October 1967 value until, for example, June 1970, 
i.e., T = 69 months. For values of T > 69, presumably the flux will 
decrease, and to a first approximation it may be assumed that this 
decrease can be estimated using the values of Rx presented in this 
section. 

As a consequence of assuming that the solar cycle parameter Rj is 
independent of B, values of Rj presented in this section may be used 
with both omnidirectional and unidirectional f luxes . 

B. The Effects of Starfish Electrons 


At some general point in time later than July 1962, the time of 
Starfish injection, the total inner belt electron flux j on a quiet day 
is 


^q ^st (18) 

where the suffixes q and st denote, respectively, the natural quiet day 
background component and the residual Starfish component. If the decay 
of the Starfish component is assumed to be exponential with character- 
istic time t, and the time dependence of the quiet day flux is written 
as 

jq " J oq £(t) (19) 
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and then the change of the total population with time is 



i.e., the difference of the source and loss terms. For epochs such that 
j » j oq > equation 20 reverts to the normal exponential decay expression. 
In the general case, however, the term j 0 q must be considered to avoid 
determining values of x too large or too small for positive and negative 
df/dt, respectively. Note that x in equation 20 represents the observed 
net depletion rate of residual Starfish electrons as opposed to the life- 
time of these particles; it is well known that cross L and pitch angle 
diffusion does occur. Hie variable x is referred to as the decay time. 

As noted in Section 6A, the major problem in determining the tem- 
poral variations is the separation of the two components, i.e., the 
residual Starfish flux and the solar cycle dependent quiet day flux. 

These components are determined in an iterative manner. Using data as 
near to the Starfish injection epoch as possible (in the present case 
September 1964), an initial decay time is determined assuming j»j 0 q 
in equation 20, and the flux is decayed using this value of T until the 
observed flux is significantly larger than the decayed flux. The dif- 
ference of these fluxes approximates to the quiet day background at this 
epoch. In practice, the epoch for this comparison is determined by the 
availability of data. With this quiet day background flux value and 
other similar values from later epochs for which it is clear that the 
Starfish residual has become insignificant, an extrapolated quiet day 
flux estimate can be made for the epochs on which the initial decay time 
is based. A corrected decay time can then be determined and the process 
repeated. The iteration procedure is lengthy and many iterations are not 
warranted, since significant error is introduced by extrapolating for 
the background flux. Further, for cases where the Starfish residual is 
small compared with the quiet day flux for all the data epochs, only a 
very approximate value of x can be determined. 
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A second iteration procedure is adopted for determining the decay 
times for the residual Starfish fluxes. In general, decay times will be 
quoted for equatorial fluxes on the assumption that T is B independent. 

It is known that this is an approximation for low L value (Stassinopoulos 
and Verzariu, 1971), although in general decay times are not determined 
for L values at which this introduces appreciable error. The pitch angle 
dependence of the residual Starfish is removed in the following manner. 

An approximate determination of the pitch angle dependence is made using 
the month of data with the greatest pitch angle range at any given L 
value. Data from all epochs are normalized to the equator assuming that 
the pitch angle distribution of the residual electrons is not time depen- 
dent, and initial values of t are determined. With these values of T, 
the time dependence is removed and, assuming that t is B independent, 
an average pitch angle distribution is obtained for all epochs. After 
normalizing to the equator, an improved T value is obtained. As with 
the previous iteration procedure, the accuracy of the technique and the 
data standard deviation do not warrant iterating more than once. In the 
following paragraphs only the final iterated values of x are presented. 

The procedure described above is adopted for analysis of the 0G0 1 
and 3 data only. No decay times are presented here for the 1963 -38C 
data, but the OGO-based values are compared with the decay time model 
of Stassinopoulos and Verzariu (1971) based on the 1963-38C data. A 
typical situation that results from separation of the two flux components 
is shown in Figure 95 for L = 1.5 and 1.9 and for 292 < E/(keV) < 690. 

The two flux components and the total flux are indicated together with 
residual Starfish flux data from OGO 1. At the lower L value, the 
dominant component at the early epochs is clearly the Starfish residual 
and, even at epoch 6/66, 51% of the total flux is Starfish residual. At 
the higher L value, the quiet day background component is discernible in 
late 1964 and early 1965, and at epoch 4/65 this component contributes 
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approximately 20% of the total flux. As a result, the final value of x 
is significantly different from the initial value obtained assuming that 
j»Joq equation 20. An approximate expression can be given for the 
ratio of the initial decay time Tf and the final value Tf accounting for 
the background component as 



where At is the time period over which xf and Tf are determined and jq 
and j st a *e evaluated at the end of the period At. For the case of 
L = 2.2, equation 21 gives that Xf = 0.81 Ti. Postulating that the 
quiet day background jq at time t may be in error by a factor of 2, the 
final decay Tf has an error from this source of ±20%. Although x f will 
be a better estimate of the decay time than tf, other error considera- 
tions such as data standard deviation dictate that the iteration pro- 
cedure is not worthwhile for j q /j s t >0.5 approximately. 

The determination of the OGO-based decay times is discussed sepa- 
rately for E < 690 keV and E > 690 keV, since in the first case the OGO 
data itself can be used for estimating the quiet day flux, whereas in 
the second case the OV3-3 equatorial perpendicular flux measurements 
must be used. The iteration procedures described in the early para- 
graphs of this section are used for E < 690 keV only. Three variables 
will be presented: the decay times, the residual flux for epoch 

September 1964, and the residual flux pitch angle dependence. 

Analysis of the low-energy data from OGO 1 indicated no residual 
Starfish fluxes except at low L values in late 1964 for 36 < E/(keV) < 
133. This is indicated in Figure 96. The broken lines on this plot 
are intended to indicate trends only. For the purposes of this diagram, 
the OGO data have been normalized to the equator using the quiet day 
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model pitch angle distribution and have been averaged. For L > 1.3, 
the general trend is for the equatorial perpendicular flux to increase 
with time. For L < 1.3, however, fluxes at the late 1964 epoch are 
above those at the later epochs by approximately a factor of 2, thus 
indicating the possible presence of a Starfish residual in 1964. The 
quality and time span of the 0G0 data at this L value are not such that 
a decay time can be determined, and it is merely noted that the Starfish 
contribution to the total flux has become insignificant by mid-1966 at 
the latest. For all other L values, all epochs covered by the OGO data 
can be used for estimating the solar cycle variation as noted in Section 
6A. 


The iterated residual flux pitch angle dependences of the Starfish- 
dominated fluxes are presented in Figures 97 through 105 for 292 > 

E/(keV) > 133 and Figures 106 through 114 for 690 > E/(keV) > 292. In 
each case the OGO 1 data have been approximately normalized to October 
1964. Simple polynomial fits to the data were made using the experssions 

logio (j s t) = a 0 + a l CQ-O q) + a2(0-ao) 2 + a3(0-ao) 3 
for otQ < G 

and (22) 


logiQ (jst) = a 0 for do > 0 

where 0 is an L-dependent cutoff for the pitch angle function analogous 
to the variable <f) used in the quiet day model (Section 2) . The L-depen- 
dent constants ao, ai, a2, a3, and 0 are presented in Table 7. By com- 
parison with Figures 14 through 57, it can be seen that the spread of the 
Starfish data is comparable to that of the quiet day data. The polyno- 
mial fits represent least squares fits weighted with the logarithm of 
the data and the standard deviation of the individual points. On each 
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of Figures 97 through 114, the Starfish flux pitch angle distributions 
are compared with the quiet day pitch angle distributions given by 
equations 6 and 7 of Section 2. In each case the quiet day distribu- 
tion is normalized at Oq = 0. At low L values , the Starfish distribu- 
tions appear flatter than the quiet day distributions. At high L 
values, a similar effect is observed, particularly if the pitch angle 
dependences are normalized at Oq = 90°, although at these L values there 
is some doubt about the value of <J>, and any differences in pitch angle 
dependence may be artificial. At intermediate L values, the two pitch 
angle dependences appear very similar in terms of both the cutoff angle 
and the variation of flux with do. Such differences as exist are ob- 
served at low ctQ, and these may be artificial because the quiet day 
pitch angle dependence is primarily influenced by the 1963-38C data in 
this region, and it was noted in Section 3 that these data indicated 
a faster rate of fall off in the cutoff region than did the 0G0 data. 

The temporal variation of the 0G0 data corrected for background and 
normalized to the equator using these pitch angle dependences is presented 
in Figures 115 to 118, Approximate decay times are given for these data 
assuming an exponential decay process. The residual fluxes are extrapo- 
lated to an epoch of June 1966 for later comparison with the quiet day 
flux at the same epoch. No decay times are given for L = 1.3, primarily 
because the data coverage is poor. In addition, however, Stassinopoulos 
and Verzariu (1971) have shown that x is B dependent in this region, and 
thus normalizing to the equator produces erroneous results. At low L 
values it is apparent that the data coverage is poor, and thus the ac- 
curacy of the decay times is suspect. Further reference is made to the 
accuracy of x later in this section. In Figure 118 the effects of the 
magnetic storm of February 7 , 1965 (maximum | Dst j = 59y) , are evident 
for L > 2.0. In addition, the equatorial perpendicular fluxes obtained 
from the 0G0 data are shown following the magnetic storm of April 18, 

1965 (maximum f Dst | = 185y). Consistent with the discussions of mag- 
netic storms presented in Section 4, no effects of the earlier storm 
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are observed for L < 2.0. No data were available for the later storm 
for L < 2.0. The injection of storm electrons at high L values evidently 
casts severe doubt on interpreting the observed flux levels at these 
energies in October 1964 and later as purely Starfish residual fluxes, 
and thus the flux is not extrapolated to June 1966 using the observed 
decay times. The variation of t from the data before and after the 
February 1965 storm is well within the standard deviation of t, and 
since storm time electrons and Starfish electrons decay in approximately 
the same manner, the observed T values are interpreted as Starfish decay 
times. A summary of the OGO-based decay times is presented in Table 8 
for 690 > E/(keV) > 292 and 292 > E/(keV) > 133. 

For the two highest energy ranges of the OGO spectrometer, little 
information is available concerning the quiet day background, and a dif- 
ferent approach to determining the decay times is adopted. In Section 
4 it was inferred from the radial profiles obtained from the high-energy 
0V3-3 data that a significant Starfish residual was present in late 1967 
at L < 1.6 approximately. A crude analysis of the high-energy OGO 1 
data is presented here in order to make approximate estimates of this 
residual flux and thus of the background flux after this component has 
decayed away. 

The pitch angle dependence of the Starfish electrons in the two 
energy ranges 1970 > E/(keV) > 690 and 4740 > E/(keV) > 1970 are shown, 
respectively, in Figures 119 through 123 and 124 through 127. These 
values were obtained by the same iterative method described above, ex- 
cept that the assumption was made that the quiet day background flux was 
small in comparison to the Starfish residual flux. Polynomial fits were 
obtained using equation 22 , and the coefficients are presented in Table 
7. In Section 4 it was assumed that the pitch angle dependence for 
electrons in the range 690 > E/ (keV) > 292 was applicable to higher 
energies. This pitch angle dependence is compared with the high-energy 
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Starfish residual pitch angle dependence in Figures 119 through 127. 

As with the lower energies, the Starfish distribution appears flatter 
at L = 1.3 and 1.4, but good agreement is observed at intermediate L 
values. This lends credence to the assumption of pitch angle distribu- 
tion energy independence for E > 292 keV. Using these pitch angle dis- 
tributions, the 0G0 1 data were normalized to the equator and are shown 
in Figures 128 and 129 for 1970 > E/(keV) > 690 and 4740 > E/(keV) > 
1970, respectively, plotted as undashed points . A limited amount of 
data was available at these energies from the 0G0 spectrometer over the 
period June through August 1966. It will be shown later that there is 
reasonable agreement between these data and the quiet day OV3-3 data 
from the same period. The averaged 0G0 and 0V3-3 data at epoch August 
1966 are shown in Figures 128 and 129, with the 0V3-3 data for epoch 
April 1967. Note that some error is involved in estimating the flux at 
energies 4740 > E/(keV) > 1970 using the 0V3-3 data because of the 
energy coverage of this satellite. In Figure 128 it is apparent that 
the Starfish component has become insignificant by August 1966 at 
L = 1.7 because solar cycle effects are evident. A similar effect is 
observed at L = 1.6 for 4740 > E/ (keV) > 1970. For higher L values, 
it will be shown that the OGO 1 data from the late 1964 period are in 
approximate agreement with the 0V3-3 data from later epochs; that is, 
the Starfish component has become insignificant by approximately early 
1965. On the assumption that the decay process is exponential, for L 
values having a significant Starfish component in April 1967 it is 
possible to make a crude estimate of the background flux and the decay 
time by removing flux increments from the total flux until the residual 
data lie in a straight line on a logarithmic plot. The results of this 
process are shown in Figures 128 and 129, where the residual flux data 
are indicated by dashed points. At all except the highest L value at 
each energy, the flux at epoch late 1964 is dominated by the Starfish 
component. At these highest L values, however, the background com- 
ponent is significant in late 1964, and the process becomes less 
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accurate because no information concerning the variation of the back- 
ground flux with solar cycle is available. It is appreciated that this 
process is highly inaccurate, and the resulting background flux values 
must be regarded as speculative. A summary of these background fluxes 
is presented in Table 9, and the decay times are shown in Table 8. 

Note that no information is presented for L = 1.3 because of the data 
coverage at this L value. In Section 4 it was shown that the magnetic 
storm effects did not affect the low L values to a significant extent, 
and these background flux values can be used for estimating the omni- 
directional flux for epochs for which the Starfish component has decayed 
to insignificance. These estimates are presented in Figure 130, and it 
can be seen that the double peak in the 500-keV radial profile is no 
longer present and that the peak flux at higher energies has moved to 
an L value of approximately 1.6. These estimated radial profiles are 
very approximate and are not incorporated into the model AE-5. 

It is interesting to compare the radial profiles of the Starfish 
and quiet day components at various epochs. These radial profiles are 
shown for epochs September 1964, June 1966, and October 1967 in Figures 
131 through 134 for the four higher energy ranges covered by the 0G0 
spectrometer. In Figure 131 for 292 > E/(keV) > 133 it is shown that 
the Starfish component has become small in comparison to the quiet day 
background by June 1966. For 690 > E/(keV) > 292, however, the two 
components are comparable in June 1966 for L < 1.6, and at L = 1.9 the 
Starfish component still contributes 20% of the total flux. By October 
1967, however, the Starfish component is clearly insignificant. A 
summary table giving the percentage contribution of the Starfish com- 
ponent to the total flux at the equator is presented in Table 10 for 
the three epochs September 1964, June 1966, and October 1967. Owing to 
the differences in the quiet day and Starfish flux pitch angle distribu- 
tions, the Starfish component generally contributes a higher percentage 
of the total flux than is indicated in Table 10 for the actual pitch 
angles used for determining the residual Starfish flux. The residual 
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flux radial profiles in Figures 131 and 132 exhibit peaks in the L range 
1.5-1. 7 depending on energy and epoch. These radial profiles are ap- 
parently not what would be expected from the artificial injection of a 
large electron population at a low L value and subsequent radial dif- 
fusion, but rather they are what would be expected from an external 
source and inward radial diffusion. The explanation lies in the varia- 
tion of the decay time with L value. It can be seen that T exhibits 
a maximum in the region of L = 1.6 -1.8 for the lower two energy ranges, 
thus resulting in a slower depletion rate at these L values. Extrapo- 
lation of these radial profiles back to epochs nearer the Starfish 
injection in July 1962 results in profiles decreasing essentially 
monotonically from low L values, as shown in Figure 135 presented for 
epoch September 1963. 

It is possible to obtain confirmation of the 0G0 residual flux 
levels by using the 1963-38C data and the OV3-3 data. Figures 136 
and 137 compare the residual Starfish flux integral spectra given by 
the 0G0 and 1963-38C satellites for epochs January 1965 and June 1966. 
The comparisons are made at B values appropriate to the 1963-38C meas- 
urements. For the data presented in Figure 136, 0G0 measurements taken 
approximately at the mirror B value were used, and thus no extrapolation 
using the polynomial fits was required. For the later epoch, however, 
some extrapolation was required. Intensity thresholds are shown for 
the 0G0 data corresponding to 1 cps. It is clear that the agreement 
between the two data sets is good. The model lines on Figures 136 and 
137 represent the integral flux predictions of the polynomial fits to 
the 0G0 data (Table 7) and the decay times t (Table 8). Some modifica- 
tions have been made to the values of ag in Table 7 to obtain the best 
fit to the two data sets shown in Figures 136 and 137. The values of 
a q in Table 7 represent the equatorial perpendicular flux with units of 
electrons/cm 2 -ster-keV and do not correspond to the arbitrary flux 
curves presented in Figures 97 through 114 and 119 through 127. Note 
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that the coefficients of the polynomials given in Table 7 should not be 
used with the decay times for decaying the flux in the energy range 
292 < E/(keV) < 690 for L > 2.0 owing to the effects of magnetic storms. 
In Figure 138, the OGO 3 and 0V3-3 data are compared for epoch August 
1966. Total integral flux is compared for those L values for which a 
Starfish residual remains in August 1966. It can be seen that the two 
data sets agree reasonably well. The averaged total flux shown in 
Figures 129 and 130 for epoch August 1966 was determined from the lines 
presented in Figure 138, and it is clear that no substantial error re- 
sults from this averaging process. For L > 1.6 the Starfish flux com- 
ponent had generally become insignificant by late 1965, as noted 
previously. The quiet day 0V3-3 and OGO data are compared for these 
L values in Figures 139 and 140 for a variety of epochs . The OGO data 
presented in these figures have been normalized to the equator using 
polynomial fits to the pitch angle distribution as with the lower L 
values. These distributions are not shown here; however, the coeffi- 
cients of the polynomial are presented in Table 7. In Figure 139 it 
can be seen that a Starfish residual remains in late 1964 at L = 1.8 
for 1970 > E/(keV) > 690. In addition, the effects of the storm of 
February 7, 1965, are evident in the OGO 1 data, and no attempt is 
made to calculate a decay time for L=1.8. At L = 1.9 and 2.0, there 
is essentially no Starfish component at these energies later than 
December 1964. Note that the agreement between the OGO 3 and OV3-3 
data in mid-1966 is not as good as at the higher L values . The OGO 
data for this epoch are based on a few measurements taken at pitch 
angles providing low count rates of order unity, and thus the 0V3-3 
measurements must be preferred. At the highest energy range of the 
OGO spectrometer shown in Figure 140, OGO data are available only for 
L = 1.7 in late 1964 and early 1965. It can be seen that a small 
Starfish residual is present in early 1965, but this effect has ap- 
parently become insignificant by late 1965. 
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From the above analysis, the times at which the Starfish component 
has become insignificant in relation to the remaining components can be 
derived. Some discussion of the word 'insignificant' is required in 
the present context. A cutoff time t c can be given by the expression 


t c = T lo &e 


Cl-P) istCt=0) 
P jq(t c ) 


(23) 


where p is the ratio of the residual Starfish to total flux at time 
t = t c . Since jq is time dependent for E < 690 keV and for T > 6, 
equation 23 must be solved by iteration. Equation 23 presupposes that 
the only significant component is the quiet day background. It has 
been demonstrated that magnetic storm effects are apparent for high L 
values and intermediate energies, and equation 23 is not applicable in 
this context since it is not possible to separate the Starfish and 
storm components. The accuracy of the foregoing analysis is such that 
solution of equation 23 is not warranted for p < 0,25 approximately. 
Equation 23 is solved for p = 0.5 (j st = jq) and p = 0.25 (jst = jq/3) 
for those L values for which storm effects are not a factor. The re- 
sulting cutoff times are presented in Table 11. The time t c is quoted 
as year and month and also as months, Tst, from Starfish injection in 
July 1962. No cutoff times are given for 133 > E/(keV) >36. It is 
clear from Figure 96 that for L > 1.3 for these energies, no significant 
Starfish residual remained later than September 1964 (T st = 26). For 
L = 1.3, an approximate cutoff date of June 1965 (T s t = 35) may be given. 
It is clear from Table 11 that the cutoff times t c show considerable 
variation with energy and L value at high energies . The longest lived 
electrons occur in the range 1970 > E/(keV) > 690 at low L values, with 
the higher energy electrons decaying away only slightly earlier. It 
should be noted that the background flux for these regions has been 
determined in a very crude manner, and errors of 6 months may be asso- 
ciated with tc. For other L values and energies, t c may be associated 
with an error of 3 months. By extrapolation from low L values for the 
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two higher energy channels, it is likely that residual electrons at these 
energies at L = 1.3 persist for longer periods. A crude estimate of t c 
at L = 1.3 can be made by extrapolating on the background fluxes given 
in Figures 133 and 134, using the averaged 0V3-3 and 0G0 3 data for epoch 
mid-1966, and by taking the decay times of Stassinopoulos and Verzariu 
(1971). Cutoff times for 690 < E/(keV) < 1970 of August 1968 and Septem- 
ber 1969 are obtained for p = 0.5 and 0.25, respectively. The corre- 
sponding dates for 1970 £ E/(keV) £ 4740 are April 1968 and March 1969. 

Little quantitative information has been given so far for the 
accuracy of the parameters determined in this section. Perhaps the most 
important variable from a practical viewpoint is the cutoff time t c . 

The accuracy of this variable depends on the value of t c and the accuracy 
of (a) the decay times t, (b) the Starfish flux model for epoch September 
1964, and (c) the solar cycle dependence of the quiet day model. While 
the accuracy of each of these three variables is not entirely independent, 
estimates of their effect on t c can be made. Taking as an example the 
case L = 1.5 and 690 > E/(keV) > 292 shown in Figure 95, if the Starfish 
model is associated with a 25% error, i.e., j s t ~ 6.92 x 10 3 ± 1.73 x 10 3 
in September 1964, this t c has an error of ±1 month for p = 0.5 and ±2 
for p = 0.25. A 25% error in the decay time is more significant in its 
effect on t c . If t = 230 ± 57 days, the error associated with t c is ±2 
months for p = 0.5 and ±4 months for p = 0.25. A 25% error in the quiet 
day flux has approximately the same effect as a corresponding error in 
the Starfish flux. The agreement between the 0G0 1 data and the 1963-38C 
data has been demonstrated (Figure 136) for values of B off the equator, 
and the pitch angle distribution is generally well defined by the OGO 
data for all oq except at high L, where a 25% error in the equatorial 
perpendicular flux is possible. However, from the above example it is 
clear that error on the decay time is the major parameter as far as the 
accuracy of t c is concerned. It is considered that the error contribu- 
tions from (b) and (c) above are relatively unimportant, except perhaps 
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at low L values and high energies, where the background flux estimate 
is associated with an error that may be as high as a factor of 2, thus 
resulting in a 6 -month error in t c . Error in the Starfish decay time r 
is a function of (1) the data time sample, (2) the standard deviation of 
the data, and (3) the significance of the quiet day background. Note 
that t is not dependent on normalizing to the equator, since it is in- 
dependent of B for the L values covered here (Stassinopoulos and 
Verzariu, 1971). As an example, if a 6-month sample of data is avail- 
able and the residual flux at either end of that period is in error by 
±10%, the resulting error in r is ±20% approximately. The decay times 
x were determined from the OGO data by weighting on the number of data 
points resulting in the monthly averaged residual flux and the standard 
deviation of the averaged figure. A representative standard deviation 
is 10% (except for the high energies where there is greater uncertainty 
in the background flux); however, the sampling period varies consider- 
ably, and thus the error on the decay time varies considerably. With 
reference to Figures 115 through 118 and 128 and 129, it can be seen 
that the sampling period varies from 4 months at low L and energy to 
2 years at low L and high energy. Estimates were made of the accuracy 
of T and are presented in Table 8. Although presented as such in this 
table, the errors on t are generally not symmetrical, and t is likely 
to be in error on the larger as opposed to the smaller side. It is 
clear that the maximum error occurs at low L values and low energies 
primarily because of the poor time coverage in this region, and these 
decay times must be interpreted cautiously. From these decay time 
errors, estimates can be made of the accuracy of the cutoff times, 
including the effects of factors (b) and (c) described at the beginning 
of the paragraph. These error estimates are given in Table 11 for 
p = 0.5 only. Because of the increased extrapolation time for p = 0.25, 
these cutoff times must be regarded as crude estimates only. 
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Verification of the OGO-based decay times can be obtained by com- 
parison with those values based on other data sets and also by compari- 
son with the decay time model of Stassinopoulos and Verzariu (1971). As 
noted previously, the decay time model is based on data from the 1963-38C 
satellite taken at Ep = 255 keV, 1.34 MeV, and 2.4 MeV. The model is 
analytic, provides T (B^Ep), and is valid in the energy range 200 < 
Ej/(keV) < 3000 approximately, t is found to be independent of B for 
L > 1.4, as noted previously. For comparison with this model, the OGO 
residual fluxes are summed to provide integral residual fluxes for 
E T = 133, 292, 690, and 1970 keV. Integral flux decay times are deter- 
mined from these residual fluxes and are presented in Table 12, together 
with errors and the corresponding values given by the Stassinopoulos 
and Verzariu model; the accuracy of the latter decay times is given as 
±15 days. Two decay times are given at L = 1.3 for the Stassinopoulos 
and Verzariu model; the high value corresponds to B = Bq, and the lower 
value to B = .22 gauss (atmospheric cutoff B c = .232 gauss approximately). 
No comparison is presented to E-p = 133 because of the energy range of 
the Stassinopoulos and Verzariu model. For the higher two channels, 
the present values are generally less than the values given by the 
Stassinopoulos and Verzariu model, although the agreement is considered 
to be reasonable. For Ep = 292 keV, however, there is substantial dis- 
agreement between the two sets of values for L < 1.8. Two legitimate 
reasons exist for differences between the two sets of decay times : (a) 
the effect of the background flux not included in the Stassinopoulos 
and Verzariu model, and (b) the variation of decay time with time. 

The present integral flux decay times are compared with those based 
on the 500 -keV data from Explorer 26 (Section 3) and those given by 
Rosen and Sanders (1971) based on the Pegasus A and B satellites for 
300 and 500 keV. It has been noted that the absolute flux levels from 
the Explorer 26 data are apparently anomalous , and decay times are de- 
termined from these data only for L values for which the quiet day 
background is clearly insignificant. The comparison between the OGO, 
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Pegasus, and Explorer 26 decay times is made in Figure 141. The 
Stassinopoulos and Verzariu model is shown as broken lines in this 
figure. The full lines are added merely to indicate trends, and the 
figures at the side of the points indicate the standard deviation of 
x. No clear pattern emerges from the comparison of the various data. 

At some energies the 0G0 decay times appear too low, at others too 
high. However, it is clear that all the data lie below the correspond- 
ing model curves by 70-140 days for Ep < 500 keV. From the comparison 
of the OGO decay times, the Stassinopoulos and Verzariu model, and data 
from other sources, a conflicting picture arises. The model indicates 
that the OGO decay times are too low by approximately 140 days for 
L < 1.8 and E-p = 292, whereas the other data indicate that the OGO 
decay times are approximately correct within the standard deviations 
of t. It is recommended that in general the most appropriate value of 
x might be the upper limit of the OGO-based values for L < 1.9 and the 
Stassinopoulos and Verzariu figures for higher L values and those regions 
where no OGO values are quoted. As noted, there is essential agreement 
between those two for Ep > 690 keV. Correspondingly, in Table 11 it is 
recommended that the upper limit of the cutoff times for channels 2 and 
3 be used. For the higher energy channels, t c in Table 11 is considered 
appropriate. The poor comparison between the various x values for 
Ep = 292 keV does cast doubt on the accuracy of the solar cycle param- 
eters for 690 > E/keV > 292 for L < 1.7 approximately, since these 
values were determined by iteration on the residual Starfish flux com- 
ponent at epoch T = 22 . Since the model decay times are greater than 
the OGO-based values , the Starfish component may be larger than that 
assumed at T = 22, and thus the solar cycle parameters presented in 
Figure 99 must be regarded as minimum values . 

As a final stage in the analysis of the effects of the Starfish 
injection, integral flux cutoff times can be determined. These are 
presented in Table 13 for p = 0.5. For comparison, three sets of decay 
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times are used: (a) the OGO-based values, (b) the values recommended 

in the previous paragraph, and (c) the Stassinopoulos and Verzariu model 
values. Points (a) and (c) can be regarded as limits on the recommended 
values given by (b). 

Earlier it was noted that equation 23 could not be used with the 
0G0 data for estimating cutoff times for those L values and energies 
for which storm effects are significant, i.e., at high L values for 
690 > E/(keV) £ 292. Integral flux cutoff times, however, can be de- 
termined for E > 255 keV using the 1963- 38C data at epochs early enough 
for the storm contribution to be small in comparison to the Starfish 
flux. Flux levels were obtained from the 1963-38C data for epoch late 
1963 and used with the decay times given by the Stassinopoulos and 
Verzariu model to determine cutoff times for p = 0.5. These are shown 
in Table 14 for 2,2 > L > 1.8. For the low L values of this range, cut- 
off times have been computed from the 0G0 data for E > 292 keV using the 
previously described method, and these are compared to the 1963-38C 
based values for E > 255 keV in Table 14. It should be noted that at 
the higher L values cutoff times have been determined using the quiet 
day background and, in relation to the background flux storm component, 
the Starfish component reaches this level of significance at earlier 
epochs than those given in Table 14. 
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7. MODEL PRESENTATION AND DATA PRODUCTS 


In this section a description of the model forms available to a 
user is given. These model forms include graphical presentation and 
a variety of computer programs that can be supplied to a user. The 
user is reminded that AE-5 is presented for an epoch of October 1967 
and that temporal variations may result in significant flux changes 
in certain regions of B-L-E space. These temporal variations include 
magnetic storm effects (Section 4) , solar cycle effect (Section 6A) , 
and residual Starfish electrons (Section 6B). The present report 
attempts to account for the last two effects. The results of these 
analyses, however, are not incorporated into computer programs and 
are presented primarily for guidance to the user. The temporal vari- 
ations, however, generally are such that the flux decreases from that 
given by AE-5 at epoch October 1967, and the model as such provides 
an adequate tool for assessing the influence that trapped inner zone 
electrons have on orbital vehicles . 

In previous documentation on trapped particle models, the major 
source of data has been in the form of omnidirectional integral flux 
tables. The present model, however, is presented in the form of 
two-dimensional carpet plots J = J(B,L) for given energy thresholds. 

A short description of the use of carpet plots is given in Appendix 
A. The omnidirectional flux is shown as a function of B and L in 
Figures 142 through 148 for E-p = 0.04, 0.1, 0.25, 0.5, 1.0, 2.0, and 
4.0 MeV. In addition, the model is presented in the more familiar 
forms of B-L flux maps and R-A plots in Figures 149 through 154 for 
energy thresholds of 40 keV, 500 keV, and 1 MeV. These three energies 
typify the three types of radial profiles in the inner belt: the 
lowest energy having a peak at L = 1.9 approximately, and the doubly 
peaked intermediate energy and the highest energy having a peak at 
L = 1.4 owing to the Starfish residual. Three-dimensional plots of 


the model flux are shown in Figures 155 through 160 for = 40, 100, 

500 keV, 1, 2, and 4 MeV. These plots provide a complete pictorial view 
of the inner belt model AE-5. The model AE-5 is incorporated into an 
orbit integration program (ORP) with the outer zone model AE-4 (Singley 
and Vette, 1972) to estimate flux doses along arbitrary orbits. The 
integration is performed for 0°, 30°, 60°, and 90® inclination circular 
orbits at altitudes in the range 150 to 18,000 n.m. These integrations 
are not presented in the present document but were reported by Singley 
and Vette (1972) with plots of flux/day as a function of altitude for 
the four inclinations. 

A wide variety of computer programs associated with AE-5 are 
available to the user. A modification of a matrix storage technique 
developed by Kluge and Lenhart (1971) for trapped radiation models 
is used in program MODEL. The program outputs omnidirectional inte- 
gral flux for an arbitrary set of input B, L and Ex values and can 
provide differential flux for a user -specified energy range. The 
outer zone model AE-4 (Singley and Vette, 1972) and a smoothed version 
of the proton models AP-5, AP-6, and AP-7 (Kluge and Lenhart, 1971; 

King, 1967; Lavine and Vette, 1969; Lavine and Vette, 1970) are in- 
cluded in program MODEL. The program may be operated with all or any 
of these models included. Versions of MODEL can be supplied for opera- 
tion on IBM 360 series, UNIVAC 1108, and CDC 6000 series computers. 

MODEL also can be supplied with plot routines suitable for use with a 
CalComp 570 plotter. The inner zone electron model provided by pro- 
gram MODEL has been smoothed for equatorial grid point integral and 
differential fluxes. This smoothing process resulted in a maximum 
5% change in the model fluxes described in the previous sections. 

The orbital integration program ORP is available to the user and 
is capable of performing integration of the model flux along an arbi- 
trary inputted orbit. Program ORP is a replacement for TRECO (Lucero, 


60 



1968) previously supplied by the National Space Science Data Center and 
requires input in the form of B, L coordinates. Programs generating B, 
L coordinates from orbital elements can be supplied by NSSDC (Kluge, 
1970; King, 1971). User documentation and program descriptions are 
supplied with MODEL and ORP (Teague et al., 1972). Model AE-5 can be 
supplied in unidirectional differential flux form for E < 690 keV and 
in integral form for 4.5 > Ej/MeV > 0.04 by program UNI. The unidirec- 
tional integral flux model can also be provided in the matrix format 
used for the omnidirectional flux model. Hie program UMODEL outputs 
unidirectional fluxes for AE-5 only. An APL (A Programming Language) 
version of the unidirectional flux model AMODEL has been included in 
the Goddard Space Flight Center library on the IBM 360/95 and can be 
used with the APL program FLUXCO to determine omnidirectional fluxes. 
Complete user documentation for programs UNI, UMODEL, AMODEL, and 
FLUXCO can be supplied (Teague, 1972). Microfilmed omnidirectional 
flux carpet plots can be supplied for user-specified energy thresholds. 

To enable the user to assess the reliability of the model AE-5, a 
system of confidence codes is presented. In developing these codes 
a number of criteria were used: number of data sets used, data coverage, 
and uncertainties introduced by temporal variations . A scale of 1 to 
10 is used, where 10 corresponds to the highest reliability with an 
expected error of a factor of 2 or less and 1 corresponds to the least 
reliability with an expected error in excess of a factor of 10. In 
general, however, efforts have been made to provide pessimistic flux 
estimates where low confidence codes are given which state that it is 
more probable that the flux is lower than the quoted value than higher. 
TVo sets of codes are given - one for the omnidirectional flux at an 
epoch of October 1967 (Table 15), and one for the integral flux solar 
cycle parameters (Table 16). In each case a brief explanation for the 
code and a section reference are given. 
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APPENDIX A. USE OF CARPET PLOTS 


Carpet plots have been used in this document for the graphical 
presentation of the solar cycle ratios (Figures 80-88 and 90-93, 

Section 6 A) and the omnidirectional model flux (Figures 142 through 
148, Section 7) . The general term for these multidimensional plots is 
nomograph (Adams, 1964), but the phrase "carpet plot" is used here to 
denote the special application for the present document. In the pres- 
ent usage, an arbitrary dependent variable Y is presented as a function 
of two independent variables Xi and X 2 . In theory, an arbitrary number 
of independent variables Xi may be covered with a practical limit of 
i equal to 4. For i > 2, however, the plots become cumbersome and 
difficult to use. 

In general, the dependent variable is plotted as a continuous 
function of an independent variable at a variety of specified values 
of the other; i.e. , 

Y = f(X lf X 2 = const) (A. 1) 

and 

Y = f (X 2 ,X : = const) (A. 2) 

Conventionally, the dependent variable is plotted as abscissa and the 
independent variables as ordinates. The origin of each ordinate scale 
is sliding, although the scales themselves remain constant. These 
scales are indicated by bars on each plot presented in this report. 

The displacement of the origin for successive curves of the family 
given by equation A.l, for example, is determined by the scale of the 
variable X 2 and the value of X 2 on those successive curves. 

PRECEDING PAGE BLANK NOT FI LIvSD 
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Carpet plots are directly amenable to interpolation. In Figure 
161 , if the value of Y is required at some point X2 = a3 and Xj = bj, 
the four points (bi,a3), (t>2> a 3)> (b3> a l)* and (b3,a2) can be easily 
located along the four lines of constant Xj and X2 given in the carpet 
plots using the scale factors for the independent variables. Other 
points lying on the curves X2 = aj and X^ = b3 can be determined in a 
similar fashion for other given curves in the two families. The re- 
quired curves X2 = a3 and X^ = b3 can then be determined to a degree of 
accuracy better than that given by simple linear interpolation, and 
the required value of Y at the point (b3, 33) can be obtained from the 
abscissa scale. 

In general, the dependent variable for a grid point, (bj, ai) in 
Figure 161 for instance, can be retrieved from a carpet plot with an 
accuracy of 1 - 2 %. Where interpolation is performed in the manner 
described above, it is estimated that this accuracy will be 2 - 3 %. 

While the omnidirectional flux data have been presented tabularly in 
previous model documentation with greater resolution than can be ob- 
tained from carpet plots, error associated with determining the inde- 
pendent variable Y from the carpet plots is considered insignificant 
in comparison to the inherent error associated with obtaining the model 
(Section 7 ) . 
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1.850257 

5 . 0 0 0 0F1 

1. 08085" 

l 

1 . 080957 

9.949556 

3.113556 

4.510856 

B . 030951 

2. 585657 

5 • 5 00 0£1 

9 . 45195" 

2 

1.353057 

1.315057 

4.505856 

6.438156 

8. 354851 

3.362057 

6.0000E1 

8.45665" 

2 

1.618657 

1.643557 

6.060956 

8.568856 

8.633651 

4.155857 

6 * 5 000E1 

7.72175" 

2 

1.872857 

1.969257 

7 . 701756 

1.080157 

8 .868351 

4.942457 

6,6000^1 

7 . 59995" 

2 

1.921957 

2.033057 

8.032656 

1.125057 

8.910051 

5.096757 

9 * 0000E1 

6.34345“ 

2 

1.921957 

2.033057 

8.032656 

1.125057 

8.910051 

5.096757 


FLUX UNITS x ELECTRONS /SQCM. SEC. STER 



TA3LE 3, SATELLITE CALIBRATION CONSTANTS 



Present Values 

Original Experimenter Values 

Satellite 

Channel 

Energy Range (keV) 

C 

Energy Range (keV) 

C 

060 1/3 

1 

36-133 

3.691 

50-120 

7.9 

( 1964-54A} 

2 

1 33-292 

6.23 

120-290 

8.4 

(1966-49A) 

3 

292-690 

6.1 

290-690 

7.2 


4 

690-1970 

6.66 

690-1700 

9.5 


5 

1970-4740 

7.794 

1 700-4000 

15.4 

Explorer 26 






(1965-78A) 

1 

>520 

2.34 x 10" 

>500 

2.5 x 10" 

19G3-38C 

1 

>254 

487.7 

>280 

600 


2 

>1340 

896 

>1200 

1200 

0V3-3 

1 

2147-2472 

- 

2147-2472 

_ 


2 

1880-2200 

- 

1880-2200 



3 

1615-1925 

- 

1615-1925 

- 


4 

1329-1651 

- 

1329-1651 

. 


5 

1075-1375 

- 

1075-1375 

. 


6 

814.5-1099,5 

- 

814.5-1099.5 

- 


7 

574.5-849.5 

- 

574.5-849.5 

- 


8 

350-600 

- 

350-600 

- 


9 

225-375 

- 

225-375 

- 


C = Flux/Count Rate (cps) 
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TABLE 4. PITCH ANGLE COVERAGE 


Data 

Time 

L Range 

cx Q Range 

Range 

B o 

8 c 

a c 

0G0 3 

6/66-1/67 

1.3 

65-80 

.17-. 145 

.1418 

.232 

51.4 


6/66-4/67 

1,4 

50-85 

.194-, 115 

.1136 

.237 

44 


6/66-11/67 

1.5 

45-80 

.185-. 095 

,0923 

.245 

37.9 


6/66-11/67 

1.6 

35-75 

.231-. 082 

.0761 

.254 

33.2 


6/66-11/67 

1.7 

30-65 

.254-. 077 

.0634 

.262 

29.5 


6-7/66 

1.8 

35-65 

.162-. 065 

.0534 

.268 

26.5 


9-11/67 

1.8 

27-50 

.259-. 091 

.0534 

.268 

26.5 


6-7/66 

1.9 

32-60 

.162-. 061 

.0454 

.274 

24 


9-11/67 

1.9 

25-45 

.254-. 091 

.0454 

.274 

24 


6-7/66 

2.0 

30-60 

. 1 56- .052 

.0389 

.28 

21.9 


9-11/67 

2.0 

22-40 

.277-. 094 

.0389 

.28 

21.9 


6-7/66 

2.2 

30-70 

.117-. 033 

.0293 

.292 

18.5 


9-11/67 

2.2 

20-28 

.25-. 133 

.0293 

.292 

18.5 ' 


6-7/66 

2.4 

20-70 

,192-. 025 

.0225 

.304 

15,8 


9-11/67 

2.4 

16-25 

.296-. 126 

.0225 

.304 

15,8 

1963-38C 

late 67* 

1.25 

61-87 

.21-. 16 

.1596 

.231 

56.5 


late 67 

1.3 

53-70 

.22-. 16 

.1418 

.232 

51,4 


late 67 

1.35-1.4 

45-63 

,23-. 16 

.1267-. 1136 

.234-. 237 

47.4-44 


late 67 

1.45 

41-51 

.24-. 17 

,1022 

.241 

40.5 


late 67 

1 .5-1.55 

37-47 

.23-. 17 

.0923-. 0837 

.245-. 249 

37.9-35,5 


late 67/late 66 + 

1.6 

35-42 

.23-. 17 

.0761 

.254 

33.2 


late 67/late 66 

1.7 

31-36 

.24-. 18 

.0634 

.262 

29.5 


late 67/late 66 

1.8 

27-32 

.26-, 19 

.0534 

.268 

26.5 


late 67/late 66 

1. 9-2.0 

22-29 

.26-. 2 

.0454-. 0389 

.274-. 28 

24-21.9 


*late 1967: 
+ 1ate 1966: 


Days 285-300 and Days 335-350 
Days 315-330 


TABLE 5. QUIET DAY PERIODS 


Jane to August 1966 
October 1966 to January 1967 
April to May 1967 
September to November 1967 
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6 


ZQUR C1QLZ ILILI MZLQS. 


L 

1.40 

1.45 

1.50 

1.55 

1,60 

1.65 

1.70 

1.75 

1.80 

1.85 

1.90 

1,95 

2.00 

2.10 

2.20 

2.30 

2.40 


R 22= FLUX (7*37)/ FLUX (7=22) 


tf ( 36-133 ) 

#(133-292) 

R ( 292-690) 

#(>2 55 ) 

1.02 

1.47 

3.30 

0.00 

1.07 

1.61 

3, 65 

0.00 

1.12 

1.75 

4.40 

0.00 

1.15 

1.83 

4,88 

0.00 

1.18 

1.92 

5.37 

2.58 

1.20 

1.99 

5 .57 

2.86 

1.21 

2,05 

5 .77 

3.15 

1.22 

2.10 

5,91 

3.28 

1.22 

2.15 

6.05 

3.42 

1.2,1 

2,17 

6,35 

3,50 

TH 

CM 

• 

H 

2.20 

6.40 

3.58 

1.20 

2.22 

6,40 

3.62 

1.20 

2.25 

6.40 

3.66 

1.20 

2.17 

6.20 

0.00 

1,20 

2.10 

6.00 

0.00 

1.20 

2.05 

5.88 

0.00 

1.20 

2.00 

5.70 

0,00 


fl = 0 DENOTES NO MEASUREMENT 



TABLE 7. COEFFICIENTS OF PITCH ANGLE DEPENDENCE OF STARFISH ELECTRONS 



292-133 keV 


690 

-292 keV 



1970-690 keV 

4740-1970 keV 

L 

a 0 

aixlO 3 

dzXlO" 

a 3 xl0 5 

0 

clo 

aixlO 3 

a 2 xl0 4 

a 3 xl0 5 

0 

a Q aixlO 3 

aaxlO" 

a 3 xlO s 

0 

a<j 

aixlO 3 a a xl0 4 

a 3 xl0 5 

n 

1.3 

4.613 

-11.94 

-1.406 

-3.291 

75 

4.447 

-21.67 

-1.007 

-4.950 

74 

3.75 -41.4 

11.9 

-11.8 

72.5 

3.057 

-32.00 .7505 

-5.181 

72.5 

1.4 

4.895 

-25.68 

-1.136 

-2.656 

71 

4.529 

2.770 

-22.14 

0 

67.5 

3.819 -22.14 

-2.603 

-.5512 

70 

3.182 

-24.37 -.7157 

-3.214 

72.5 

1.5 

4.959 

.4001 

-5.455 

0 

77.5 

4.622 

-20.38 

5.598 

-5.702 

66 

3.750 -14.46 

-13.44 

0 

67 

2.857 

-14.33 -3.551 

-.6104 

74.5 

1.6 

4.989 

-7.791 

3.874 

-2.032 

74 

4.549 

-22.23 

-2.539 

-1.328 

64.5 

3.494 -9.510 

-9,798 

0 

67 

2.321 

-15.69 -5.382 

-.4008 

67 

1.7 

5.004 

-13.06 

1.852 

-1.063 

75 

4.412 

-18.90 

-5.672 

-.6142 

63.5 

3.037 -5.956 

-9.806 

0 

66 





1.8 

4.951 

-15.28 

-3.408 

-.3547 

65 

4.250 

-22.19 

-6.253 

.1122 

60.5 

2.629 -16.67 

-5.856 

-.0271 

65 





















NO DATA 



1.9 

4.913 

1.128 

-8.49 

0 

70 

4.033 

-11.82 

-2.910 

-1.070 

63.5 

2.255 -16.59 

-6.773 

.4596 

65 





2.0 

4.839 

-7.729 

1.624 

-.8415 

79 

3.792 

-18.58 

7.626 

-1.882 

75 

1.998 -16.32 

-5.651 

-.4836 

65 





2.2 

4.586 

-2.266 

-2.219 

-.4664 

70 

3.29 

-4.455 

-1.576 

-.2088 

80 

NO DATA 








Coefficients provide flux in units of electrons/cm 2 sec ster keV at epoch September 1964. 




TABLE 8. STARFISH DECAY TIMES 



Energy Ranges (keV) 

L 

292-133 

690-292 

1970-690 

4740-1970 

1.4 

163 ± 60 

216 ± 70 

380 ± 20 

310 ± 20 

1.5 

200 ± 70 

230 ± 70 

375 ± 20 

310 ± 20 

1.6 

233 ± 45 

260 ± 50 

312 ± 20 

290 ± 20 

1.7 

265 ± 42 

230 ± 42 

310 ± 40 

ND 

1.8 

267 ± 37 

230 ± 37 

ND 

ND 

1.9 

250 ± 30 

211 ± 30 

ND 

ND 

2.0 

219 ± 30 

150 ± 25 

ND 

ND 

2.2 

98 ± 30 

103 ± 25 

ND 

ND 


Decay times are in days. 
ND denotes no data. 








UUK 9 2I£Z£ KKK&Ql MQKQKQUHa ILU ££II{*dZ££ 

ENERGY RANGE ( KEV ) 


h 

690-1970 

1970-4740 

1.4 

17 9.8 

27.2 

1.5 

249.7 

35.8 

1.5 

299.7 

41.3 


EQUATORIAL PERPENDICULAR FLUX ELECTRONS /SQCM .SEC . STER.KEV 
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X42££ 10 S.LQELEIQ.AHQE. QZ £2 A&IIS.K ILUl QQdE2MUl 


STARFISH COMPONENT AS PERCENTAGE OP TOTAL FLUX 


SEPTEMBER 1964 


CR2 

CR 3 

CR 4 

CHS 

95 

99 

ND 

ND 

75 

98 

97 

98 

72 

98 

95 

96 

71 

98 

90 

86 

69 

97 

82 

65 

67 

95 

70 

0 

66 

92 

0 

0 

63 

89 

0 

0 

59 

81 

0 

0 


JULX 1966 


CN 2 

CR 3 

CRk 

CR 5 

ND 

ND 

ND 

ND 

4 

52 

87 

87 

8 

52 

75 

71 

10 

51 

56 

40 

12 

36 

29 

0 

11 

24 

0 

0 

9 

15 

0 

0 

6 

ND 

0 

0 

0 

ND 

0 

0 


OCTOBER 1967 


CR2 

CR 3 

CR* 

CHS 

ND 

ND 

ND 

ND 

0 

3 

68 

60 

0 

3 

54 

47 

0 

3 

24 

5 

0 

1 

8 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


CHANNEL 

2 

3 

4 

5 


ENERGX RANGE ( KEV ) 
133-292 
292-690 
690-1970 
1970-4740 


ND DENOTES NO DATA. 

0 DENOTES LESS THAN 1 PERCENT 



TABLE 11. STARFISH FLUX CUTOFF TIMES 



P 

= 0.5 



P = 

0.25 


L 

Ch 2 

Ch 3 

Ch 4 Ch 5 

Ch 2 

Ch 3 

Ch 4 

Ch 5 

1.4 

2/65 (31) ±3 

7/66 (48) i6 

3/68 (68) ±8 1/68 (66) +7 

8/65 (37) 

11/66 (52) 

7/69 (84) 

12/68 (77) 

1.5 

3/65 (32) ±3 

7/66 (48) ±6 

9/67 (62) +7 4/67 (57) ±6 

9/65 (38) 

11/66 (52) 

10/68 (75) 

3/68 (68) 

1.6 

3/65 (32) ±2 

7/66 (48) ±4 

9/66 (50) ±2 3/66 (44) ±2 

10/65 (39) 

11/66 (52) 

8/67 (61) 

2/67 (55) 

1.7 

4/65 (33) ±2 

2/66 (43) ±3 

2/66 (43) ±2 -10/65 (39) ±2 

12/65 (41) 

8/66 (49) 

8/66 (49) 

-4/66 (45) 

1.8 | 

3/65 (32) ±2 

2/66 (43) ±3 

ND <12/64 (< 29) 

10/65 ( 39) 

6/66 (47) 

ND 

<12/64 (< 29) 

1.9 

1/65 (30) ±2 

11/65 (40) ±3 

ND <12/64 (< 29) 

9/65 (38) 

4/66 (45) 

ND 

<12/64 (< 29) 

2.0 

1/65 (30) ±2 

ND 

ND <12/64 (< 29) 

7/65 (36) 

ND 

ND 

<12/64 (< 29) 

2.2 j 

10/64 (27) ±2 

NO 

ND <12/64 f: 29) 

2/65 (31) 

ND 

ND 

<12/64 (< 29) 


NO denotes no data because of magnetic storm effects or no measurements. 
Figures in parentheses represent months from Starfish injection. 


Channel 

2 


Energy Range (keV) 
133-292 




TABLE 12. COMPARISON OF INTEGRAL FLUX DECAY TIMES 




Energy Threshold (keV) 


133 

292 

690 

1970 

L 

Present 

Present 

S and V* 

Present 

S and V* 

Present 

S and V* 

1.3 

ND 

ND 

382-333 

ND 

368-317 

ND 

325-265 

1.4 

230 ± 60 

265 ± 70 

415 

348 ± 20 

390 

310 ± 20 

340 

1.5 

245 ± 60 

280 ± 70 

445 

350 ± 20 

410 

310 ± 20 

320 

1.6 

256 t 45 

270 ± 45 

445 

305 ± 20 

400 

290 ± 20 

270 

1.7 

240 ± 42 

240 ± 42 

420 

310 ± 40 

360 

ND 

180 

1.8 

250 ± 37 

230 ± 37 

360 

ND 

300 

ND 

80 

1.9 

240 ± 30 

211 ± 30 

290 

ND 

215 

ND 

20 

2.0 

200 ± 30 

150 ± 25 

200 

ND 

135 

ND 

ND 

2.2 

100 ± 30 

103 ± 25 

64 

ND 

15 

ND 

ND 


Decay times are in days. 

ND denotes no data. 
*Stass1nopoulos and Verzariu, 



TABLE 13. COMPARISON OF INTEGRAL FLUX CUTOFF TIMES 


Energy Threshold (keV) 



133 

292 

690 

1970 

L 

a 

b 

a 

b 

c 

a 

b 

c 

a 

b 

c 

1.4 

10/63 (39) 

1/66 (42) 

7/66 (48) 

11/66 (52) 

2/67 (55) 

4/68 (69) 

5/68 (70) 

9/68 (74) 

2/68 (67) 

4/68 (69) 

4/68 (69) 

1.5 

9/65 (38) 

11/65 (40) 

6/66 (47) 

9/66 (50) 

12/66 (53) 

8/67 (61) 

10/67 (63) 

2/68 (67) 

3/67 (56) 

4/67 (57) 

4/67 (57) 

1.6 

7/65 (36) 

8/65 (37) 

5/66 (46) 

7/66 (48) 

10/66 (51) 

8/66 (48) 

9/66 (49) 

2/67 (55) 

6/66 (44) 

4/66 (45) 

4/66 (45: 

1.7 

6/65 (35) 

8/65 (37) 

1/66 (42) 

3/66 (44) 

7/66 (48) 

9/65 (38) 

10/65 (39) 

11/65 (40) 

ND 

-10/65 (39) 

ND 

1.8 

5/65 (34) 

6/65 (35) 

8/65 (37) 

10/65 (39) 

12/65 (41) 

ND 

ND 

ND 

ND 

ND 

ND 

1.9 

3/65 (32) 

4/65 (33) 

4/65 (33) 

5/65 (34) 

6/65 (35) 

ND 

ND 

ND 

ND 

ND 

ND 

2 

2/65 (31) 

3/65 (32) 

ND 

ND 

ND 

ND 

ND 

ND 

ND 

ND 

ND 

2.2 

10/64 (27) 

10/64 (27) 

NO 

NO 

ND 

ND 

ND 

ND 

ND 

ND 

ND 


NO denotes no data. 

a. 0G0 decay times. 

b. Recommended decay times. 

c. Stassinopoulos and Verzariu model decay times. 



r42£ £ 14 


QVZ Q ll HUES. EEQM 196 3-38 C QdZd 



1963-38C 

OGO 1 

L 

ET- 255 KEV 

ET* 292 KEV 

1.8 

3/66 (44) 

12/65 (41) 

1.9 

9/65 ( 3 B ) 

6/65 (35) 

2 « 0 

3/65 (32) 

ND 

2.2 

6/64 (23) 

ND 


ND DENOTES NO DATA . 

FIGURES IN PARENTHESIS DENOTE MONTHS FROM STARFISH. 



TABLE 15. OMNIDIRECTIONAL FLUX CONFIDENCE CODES 


Code 

B Range 

L Range 

Ej Range 

Section 

Conment 

1 

> B o 

1.2-1. 4 

>3 MeV 

4 

Extrapolation on both B dependence and 
spectrum, no data 

2 

~ B o 

1.2-1* 4 

>3 MeV 

4 

Extrapolation on spectrum, no data 

4 

> B o 

1. 2-1*7 

>250 keV 

6B 

Possible presence of Starfish electrons 

5 

> B o 

1.9-1. 4 

4-2 MeV 

4 

Magnetic storm effects, single data set, 
B extrapolation 

6 

>B o 

1*7-1. 9 

>500 keV 

4 

Single data set, B extrapolation 

6 

lB o 

< T. 25 

all energies 

5 

L extrapolation 

6 

~ B o 

>1*5 

all energies 

3 

Poor data 

7 

iB o 

1.3 

all energies 

3 

Poor 0G0 data 

8 

V 

V 

CD 

O 

>2 

all energies 

3 

Poor pitch angle coverage 

10 

iB o 

1 .4-1.9 

< 250 keV 

3 

Agreement between three data sets 
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TABLE 16. INTEGRAL FLUX SOLAR CYCLE PARAMETER CONFIDENCE CODES 


Code 

L Range 

E r 

T Range 

Section 

Comment 

3 

<1.8 

250,500 

>22 

6B 

Significant Starfish flux at T»22 re- 
sulting in iteration 

4 

>1.9 

all 

all 

3 

Poor OGO data at high L values 

4 

all 

40 

all 

6A 

Small Ry values; data standard deviation 
becomes significant. 

5 

<1.4 

250,500 

all 

6B 

Hardening of spectrum; assumed constancy 
of j (>690) term in equation 16 

5 

>1.8 

500 

all 

6A 

Storm effects term in equation 17 becomes 
significant 

7 

1. 6-2.0 

250,500 

all 

6A 

Two data sets available (i.e., OGO and 
1963-38C) 


Note that these confidence codes are low because integral flux values of Rj are determined from the 
OGO data using an approximate expression, and B independence has been assumed, Further, if is 
used to extrapolate beyond the epoch of October 1967, as described in Section 7, the above confi- 
dence codes will be reduced because of asymmetries in the solar cycle. 
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(Jj ) r <cm 2 -sec-ster) - 



Figure 1, Radial Profile of the Reference Electron Content Parameter 


(Jj) r approximates to the total unidirectional electron flux above 
thermal energies at the reference pitch angle a tt = 90° (Section 2). 


PRECEDING 1 



^ ft ***** 



_ — T> 










/ / / Jensen and Cain, 1960 


/// h min = 100 km 

// / X. Model B c 



1.4 1.6 1.8 2.0 2.2 2.L 


L Value (earth radii) 
Figure 4. Atmospheric Cutoffs B, 




FIGURE 6 
QUIET DAY MODEL 
PERPENDICULAR INTEGRAL FLUX 
GREATER THAN 4 0 KB V 
EPOCH OCTOBER 1967 



FLUX units - uicTtom/BQ CM ABC am 


FIGURE 6 
QUIET DAY MODEL 
PERPENDICULAR 1NTEORAL FLUX 
GREATER THAN 100 KEV 
EPOCH OCTOBER 1967 


F 1 0 II R E 7 
QUIET DAY MODEL 
PERPENDICULAR INTEGRAL FLUX 
GREATER THAN 250 KEV 
EPOCH OCTOBER 1967 


FIGURE a 


GU1BT DAY MODEL 
ELECTRON CONTENT PARAMETER 
EPOCH OCTOBER 1 V « T 




1963-38C 


1964 


■« » 
Explorer 26 


OGO 1 
Low L Value 


OGO 1 High L 
Value 


Jan. 1965 


OV3-3 


OGO 3 Low L 
Value 


OGO 3 High L Value 


Jan. 1966 Jan. 1967 Jan. 1968 
Time 


Figure 10. Data Coverage 



2 03 


244 0.7452 200-1400 4 

:xp 254 0.7323 40-120 4; 

exp 256 0.7429 40-100 4 

Geom. Factor = 2.8 x 10' 3 cm 2 -ster 


Flux j(>E T ) = (cps) x C 


04 05 06 07 08 09 10 11 12 1 



Hardness Parameter Xj^ or X 2 (keV) 
Figure 11 . 1963-38C Calibration Curves 






Geom. Factor = 2.8 x 10 3 cm 2 -ster 
Flux j (>E t ) = (cps) x C 



4 5 6 7 8 9 1 


E 0 x 10 2 (keV) 


Figure 12. 1963-38C Calibration Curves 








Figure 13. Explorer 26 Calibration Curves 
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The radial profile given by the 0V3-3 data in August 1966 is compared 
with the analytical quiet day flux description given in Section 2 
(broken line). A residual Starfish flux is observed for L s 1.8 at 
this epoch, and the Starfish model (Section 6B) is used to estimate 
the total flux in August 1966 tsolid line). 
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Figures 16-57 . These computer-generated plots compare the present 
inner zone electron model with the 0G0 1, 0G0 3, and 1963-38C data. 

The 0G0 data are presented at a variety of epochs to demonstrate 
solar cycle effects. Each model epoch is identified by a value of 
T measured in months from solar minimum taken as September 1964. 

T = 37 corresponds to the model epoch of October 1967. The ordinate 
scale on each plot is arbitrary counts/sec, and the conversion factors 
for flux are presented in Table 3. The 0G0 data are monthly averaged 
values from the time periods indicated in the legends. A consistent 
set of symbols is used: 

A September to November 1967 (T = 37 approx.) 

□ April to May 1967 (T = 31 approx.) 

* December 1966 to January 1967 (T = 27 approx.) 

+ June 1966 to July 1966 (T = 22 approx.) 

X September 1965 (T = 12 approx.) 

0 September 1964 to March 1965 (T = 6 approx.) 

For the lower L values and for the 292- to 690-keV data, the estimated 
quiet day flux after removal of the Starfish residual flux is shown 
for T = 22 labeled as "T = 22 EST BACK." The determination of the 
solar effects and Starfish residual fluxes is discussed in Sections 
6A and 6B, respectively. On each 0G0 plot the atmospheric cutoff 
equatorial pitch angle is given in degrees (Section 2). 

Data from the 1963-38C satellite are normalized to an epoch of late 
1967 (T = 37 approx.). The legend indicates the time period of the 
actual measurements as year and day range in' that year. A consistent 
set of symbols is used: 

+ 1967 Day 335-350 

A 1967 Day 285-300 

□ 1966 Day 315-330 

For 1.3 < L/ (earth radii) <1.4, "TOTAL" and "QUIET DAY" model lines 
are shown for T = 37. The differences arise because of the hardening 
of the spectrum at low L values for E > 690 keV (Section 5A) . On each 
1963-38C plot, the atmospheric cutoff B value is given in gauss. 
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0GQ3/1 292-690 KEV, L=1.4 FIG. 24 
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ARBITRARY COUNT RRTE. FIGURES 32 - 35 


0G03/1 36-133 KEV, 1=1.6 FIG. 32 



0G03/1 133-292 KEV. L=1.6 FIG. 33 



0GO3/1 292-690 KEV. L=1.6 FIG. 34 


1963-38C 255 KEV. L=1.6 FIG. 35 
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QGQ3/1 36-133 KEV, 1=1.8 FIG. 40 



0GG3/1 133-292 KEV, L=1.8 FIG. 41 



0G03/1 292-690 KEV. L=1 . 8 FIG. 42 



1963-38C 255 KEV. L=1.8 FIG. 43 
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1963-38C 255 KEV. 1=1.9 FIG. 47 
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1963-38C 255 KEV, L=2.0 FIG.51 
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Figure 60. Quiet Day High-Energy Equatorial Perpendicular Flux 

The 0V3-3 quiet day differential spectra are shown for high L values. 
An increase in flux is observed at the lower energies from the earlier 
(broken line) to the later epoch (solid line) for L = 1.9 and 2.4. A 
similar effect (not shown) is observed at L = 2 and 2.2. 
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Figure 61. Quiet Day High-Energy Equatorial Perpendicular Flux 


The 0V3-3 quiet day differential spectra are shown for selected L val- 
ues below 1.9 for two epochs. At L = 1.3 Starfish decay is observed, 
while at L = 1.8 an increase in the flux with time is observed, thus 
indicating solar cycle effect. 
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Figure 64. Storm Model 
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Figure 66. Probability of Storm Flux Contribution 

The probability of the flux being greater than the quiet day value 
(solid line) or the averaged storm time value (broken line) is shown for 
L = 1 .9 and L = 2.8, 
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Figure 68. Modifications to Quiet Day Model 

The high-energy integral flux spectra are shown for quiet times (broken 
line) and for the averaged storm time situation (solid lines) indicat- 
ing the increased importance of storm effects at the higher L values. 
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Figure 69. Modifications to Quiet Day Model 

The high-energy integral flux spectra are shown for quiet times (broken 
line) and for the averaged storm situation (solid lines). At L = 1.7 
magnetic storm effects do not significantly influence the average model 
and no modification is made to the quiet day flux. 
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Figure 70. Modifications to Quiet Day Model 


The low-energy integral flux spectra are shown for L ^ 1.6. The modi- 
fications to the low-energy quiet day flux based on the low-energy 0G0 
data (broken line) are shown resulting from the hardening of the spec- 
tra at higher energies. 
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Figure 71. High-Energy Differential Spectra L <1.6 

Differential spectra are shown for low energies and low L values. The 
solid lines indicate the analytical quiet day flux and the broken lines 
indicate that monotonic spectra are given by the 0V3-3 data for higher 
energies. 
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Figure 72. Extrapolated High-Energy Spectra 
Based on 0V3-3 Data 
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Figure 73. Extrapolated High-Energy Spectra 
Based on 0V3-3 Data 
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Figure 86. OGO-Based Solar Cycle Ratios R^» 36 < E/(keV) < 133 
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Figure 87. OGO-Based Solar Cycle Ratios Ry, 292 ;> E/(keV) £ 133 
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Figure 88. OGO-Based Solar Cycle Ratios R T , 690 > E/(keV) a 292 
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Figure 91. Integral Flux Solar Cycle Ratios Ry, E > 100 keV 
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Figure 92. Integral Flux Solar Cycle Ratios Rj, E > 250 keV 
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Figure 94. Sunspot Number R z 

The smoothed Zurich sunspot number is shown as a function of time in- 
dicating that the solar cycle parameters have been obtained over a 
period covering most of the change from solar minimum to solar maxi- 
mum. The numbers in parentheses on the abscissa scale are months from 
solar minimum. 
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Figure 95. Separation of Flux Components 

The relative importance of the Starfish and quiet day flux components 
is shown as a function of time for the interval 690 < E/(keV) s 292 
and L = 1.5 and 1.9. Representative data only are shown for the June 
1966 and October 1967 epochs. The numbers in parentheses on the ab- 
scissa scale are months from solar minimum. 
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Figure 96. Temporal Variation of Low-Energy 060 Data 

The variation of count rate with time is shown for the lowest energy 
of the OGO spectrometer. The broken lines indicate trends, and the 
possible presence of Starfish electrons at L = 1.3 is evident. 




Figures 97-114 and 119-127 . These computer- generated plots present 
comparisons between the residual Starfish electron flux data from the 
OGO 1 satellite and the Starfish model discussed in Section 6B. The 
ordinate scale is arbitrary residual Starfish counts/sec and may be 
converted to flux using the constants given in Table 3. On each plot 
comparison is made between the polynomial fit of the Starfish model 
(labeled POLY FIT) and the quiet day model normalized to the equatorial 
flux given by polynomial fit (labeled QUIET DAY) . The OGO data and 
the model lines are normalized to an epoch of September 1964. On each 
plot the atmospheric cutoff equatorial pitch angle is given in degrees . 


145 



ARBITRARY CPS ARBITRARY CPS 


RESIDUAL STARFISH FLUX FIGURES 97-100 


QGG 1 133-292 KEV, L=1.3 FIG. 97 



OGQ 1 133-292 KEV t L=1.4 FIG. 98 



ALPHA 0 (DEGS) 


OGG 1 133-292 KEV, L=1.5 FIG. 99 OGO 1 133-292 KEV, L*1.6 FIG. 100 




146 






ARBITRARY CPS ARBITRARY CPS 


RESIDUAL STARFISH FLUX FIGURES 101 - 104 


0G0 1 133-292 KEV, L=1.7 FIG. 101 



OGQ 1 133-292 KEV, L=3.8 FIG. 102 



OGO 1 133-292 KEV f L=1.9 FIG. 103 



ALPHA 0 (OEGSJ 


OGO 1 133-292 KEV. L=2.0 FIG. 104 



147 





ARBITRARY CPS ARBITRARY CPS 


RESIDUAL STARFISH FLUX FIGURES 105 - 108 


0G0 1 133-292 KEV, L=2.2 FIG. 105 



GGQ 1 292-690 KEV. L=1 . 4 FIG. 107 
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Figure 115. Temporal Variation of Residual Starfish Flux 

Based on 0G0 1 Data 
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Figure 116. Temporal Variation of Residual Starfish Flux 

Based on OGO 1 Data 




Figure 117. Temporal Variation of Residual Starfish Flux 

Based on OGO 1 Data 
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Figure 128, Temporal Variation of High-Energy Residual Starfish Flux 
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Figure 129. Temporal Variation of High-Energy Residual Starfish Flux 




Figure 130. Estimated High-Energy Omnidirectional Background Flux 

The AE-5 model radial profiles are shown for a variety of energies 
(solid lines) for epoch 1967. The broken lines indicate the estimated 
flux after the Starfish flux component has decayed completely. 
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Figure 131. Flux Component Radial Profiles 
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Figure 133. Flux Component Radial Profiles 

The Starfish flux and quiet day flux components are compared. The 
broken line indicates the estimated flux after the Starfish flux com- 
ponent has decayed completely. No time dependence is given for this 
quiet day flux, and thus it has an arbitrary epoch. 
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Equatorial Perpendicular Count Rate 



Figure 134. Flux Component Radial Profiles 


The Starfish flux and quiet day flux components are compared. The 
broken line indicates the estimated flux after the Starfish flux com- 
ponent has decayed completely. No time dependence is given for this 
quiet day flux, and thus it has an arbitrary epoch. 
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Figure 135. Extrapolated Starfish Radial Profiles 

The Starfish model radial profiles for September 1964 are extrapolated 
to September 1963 to indicate the monotonic radial profile for this 
epoch. 
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Figure 136. Comparison of 0G0 3 and 1963-38C 
Integral Residual Spectra 
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Figure 137. Comparison of 060 3 and 1963-38C 
Integral Residual Spectra 
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Figure 138. Comparison of 0G0 3 and 0V3-3 Integral Spectra 
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Figure 140. Comparison of 0G0 3 and 0V3-3 Flux 
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Figure 141. Comparison of OGO, Pegasus, and Explorer 26 Decay Times 

Starfish decay time data from various satellites are compared with the 
Stassinopoulos and Verzariu model (broken lines). The solid lines in- 
dicate the trends suggested by the data. The figures at the side of 
the data points represent estimates of the errors in days. 



Figures 142-148 . These computer-generated plots present carpet plots 
of the AE-5 omnidirectional flux as functions of B and L for threshold 
energies Ej = 0.04, 0.1, 0.25, 0.5, 1.0, 2.0, and 4.0 MeV. A descrip- 
tion of the use of these carpet plots is given in Appendix A. In 
general, lines of constant B are presented in 0.02-gauss increments 
from the equator to 0.28 gauss, and lines of constant L are presented 
in increments of 0.05 earth radii for 1.2 < L < 1.6 and increments of 
0.1 earth radii for 1.6 < L < 2.4. In some cases, lines are omitted 
for clarity. For the energies 1.0, 2.0, and 4.0 MeV, the plots are 
subdivided into two L ranges because of steep gradients in the radial 
profiles at these energies. In each figure the ordinate scale incre- 
ments are shown as error bars on the left-hand side of the plot and the 
abscissa scale is shown as powers of ten. 
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FIGURE 142 

AE5 OMNIDIRECTIONAL INTEGRAL FLUX, ET=40 KEV 
EPOCH OCTOBER 1967 
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FIGURE 143 

AE5 OMNIDIRECTIONAL INTEGRAL FLUX, ET=tOO KEV 
EPOCH OCTOBER 1967 
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FIGURE 144 

AE5 OMNIDIRECTIONAL INTEGRAL FLUX, ET=250 KEV 
EPOCH OCTOBER 1967 
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FIGURE 145 

AE5 OMNIDIRECTIONAL INTEGRAL FLUX, ET=500 KEV 
EPOCH OCTOBER 1967 
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FIGURE 146 

AE5 OMNIDIRECTIONAL INTEGRAL FLUX, ET= 1 MEV 
EPOCH OCTOBER 1967 
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FIGURE 14 6 CONT 

AE5 OMNIDIRECTIONAL INTEGRAL FLUX, ET= 1 MEV 
EPOCH OCTOBER 1967 
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FIGURE 147 

AE5 OMNIDIRECTIONAL INTEGRAL FLUX, HT=2 MEV 
EPOCH OCTOBER 1967 
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FIGURE 147 CONT 

A 155 OMNIDIRECTIONAL INTEGRAL FLUX, ET=2 MEV 
EPOCH OCTOBER 1967 
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FIGURE 148 

AE5 OMNIDIRECTIONAL INTEGRAL FLUX, ET=4 MEV 
EPOCH OCTOBER 1967 
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FIGURE 148 CONT 

AE5 OMNIDIRECTIONAL INTEGRAL FLUX, ET=4 MEV 
EPOCH OCTOBER 1967 
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FIGURE 166 

ABB OMNI-DIRECTIONAL INTBORAL FLUX 
ORBATER THAN 100 KEV 
EPOCH OCTOBER 19BT 
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FIGURE 


A B 5 OMNI-DIRECTIONAL I HTBGRAL FLUX 
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